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Abstract 


Theoretical foundations of our understanding of ‘the problem of the 
termination of the solar wind are reexamined in the light of most recent 
findings concerning the states of the solar wind and the local interstellar 
medium. The investigation suggests that a simple extention of Parker’s 
(1961) analytical model provides a useful approximate description of the 
combined solar wind, interstellar wind plasma flowfield under conditions 
presently thought to occur. A linear perturbation solution exhibiting 
both the effects of photoionization and charge exchange is obtained for 
the supersonic solar wind. The solution demonstrates the point that the 
addition of mass and decrease of velocity as the result of photoionization 
in the supersonic region of the solar wind has no effect on the momentum 
flux and, thus, the location of the shock surface. The effect of charge 
exchange, however, serves to reduce by about 147o the estimate of the 
shock radius provided by Farker's solution as extended. Finally, investi- 
gation along lines of the one dimensional compressible flow equation and 
heating associated with polytrope solutions brings to light further theo- 
retical evidence that the classical coefficients of thermal conductivity " 
for an ionized gas may be too high - perhaps by as much as a factor of 
five to describe the flow of the supersonic solar wind beyond earth. 

Derived from its roots in kinetic theory and the plasma transport 
equations a ^ Braginskii (1965), a model for the combined flow of the 
solar wind with the ionized component of the interstellar wind ijs formu- 
lated on the level of one-fluid magnetohydrodynaraics . The mode,! features 

/ 

quasi-normal, strong shock relations on a closed surface surrounding the 
sun and imbedded in the solar wind. Tangential discontinuity relations 
apply at the free surface between the solar wind and ionized interstellar 



Chapter I. - The Problem 


1, Introduction - Physical and Historical Perspective 

From the theoretical point of view the termination of the solar wind 
is the problem of finding a model for the extended solar corona such that 
admissible boundary conditions match conditions in the perturbed local 
interstellar medium. Taken in its full physical and attendant mathemati- 
cal generality the problem is so complex that the pioneer analyst in the 
subject Parker (1963) quickly and, by events, correctly pronounced that 
”real progress can come only ^ hoc ". 

After a brief statement of conditions in the local interstellar 
medium as best they are known to astronomers, we bring forth in this 
introductory section two such ad hoc models of the corona. The models 
which combine the virtues of mathematical tractability with broad pheno- 
menological richness serve to illuminate the character of the problem 
and exhibit important qualitative and semi -quantitative features of its 
solution. The models are the static, heat conducting corona and the 
dynamic polytrope corona. 

The former model has its roots in the earliest description (Chapman, 
1957) of the far corona but a solution satisfying a boundary condition of 
greater physical content has been found by the present author and is given 
here. The latter model incorporates the prediction of the interplanetary 
supersonic solar wind (Parker, 1958). The material, which includes the 
critical polytrope solution chosen by Parker, is treated in the concise, 
explicit manner due to Bondi (1952) who had invoked the other member of 
the family of critical polytrope solutions earlier in the study of stellar 
acretion, 

Axford (1972), Wentzel (1972), and Silk (1973) have recently sum- 
marized theory and observation concerning the state of the interstellar 
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medium in the vicinity of the sun. Neglecting infrequent small dense 
clouds whose inhomogeneous effect is not the subject. of concern here, 
the picture that emerges is of a tenuous, partially ionized, weakly mag- 
netized, hot hydrogen gas. Typical numbers on which the discussion is 
based herein appear in Table 1. 


Table 1. Properties of the 

local interstellar medium. 

TEMPERATURE 

10^ - 10^ ° K 

^ARTICLE DENSITY 



NEUTRAL HYDROGEN 

0.1 - 0.2 cm 

PROTON OR ELECTRON 

0,01 - 0.025 cm 

MAGNETIC FIELD STRENGTH 

—6 

3 X 10 ” gauss 


a. Static, heat conducting corona 

*^hc model accumcc a spherically synnietric corona of fully ioni^ed- 
hydrogen gas. The equation of static equilibrium yields 

42 = - GM m (l.la-1) 

dr o 

r 

wti^ere at r the radial distance from the center of the sun p is the 
pressure, m is the mass of a proton and an electron, n is the number 
density of protons or electrons, is the solar mass, and G is the 

universal gravitational constant. 

Starting at some high temperature T^ at very close to 

a constant heat flux is assumed satisfying the energy equation 




div (q) = 0 


(l.la-2) 
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where q is the conduction heat flux density vector given by 



(l.la-3) 


In equation (3) the coefficient of thermal conductivity k appropriate 
to fully ionized hydrogen can be expressed as 


K = A T 


5/2 


(l.la-4) 


where A is a constant whose value as computed by Chapman (1954) 
is 4 X 10^^ ill CSS units when the temperature is measured in K . 

With the spherically symmetric form of the divergence operator, 
equation (2) with equations (3) and (4) provides the &rst integral 


7 S/2 
-Ar T ' 


dT 

dr 


‘o o 


<r.la«5J 


Equation (5) is integrated over the interval r^ to 


r yielding 



(l.la-6) 


Letting the temperature asymptotically approach the temperature of 

7 

the interstellar medium permits the evaluation of the constant ^ ^ 

in equation (6). (Here Chapman (1957) arbitrarily took the temperature 
to vanish at infinity,) 

As the result, equation (6) may be represented 


T 



+ 




(iaa-7) 
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Solving for T we have 


T := T 




2/7 


/T v 7/2 r > 

i+(r) (7 -1) 


2/7 


(l.la-8) 


When T « T equation (8) predicts that for a vary wide range in r 

such that — « j the temperature in the corona drops off as 

-2/7 

r 

Equation (8) and the equation of state 


p = 2 k h T 


(l.la-9) 


permit the equation of equilibrium (1) to be cast 


dr 


- e 


T 7/2 -2/7 5/7 -12/7 

7 1 '• ' 


(1. la-10) 


where € = 


GM m 

O : 

2kT r 
o o 


/T„\ 7/2 

and X = 1 - {y) 


With the change of variable x' = t/r^ > equation (10) may be rewritten 
in the form 


p 




7/2 -2/7 


dr^ 

.2 


(1. la-11) 


Equation (11) is integrated over the interval 1 to r"* yielding the 


solution 




In the limit as r -♦ « , 



The corresponding number density is 



(1. la-13) 


(1. la-14) 


Estimating r^ by = 7 X 10^^ cm the radius of the visible 

4o ^-^6o. 

disc and taking T =10 K from Table 1, ± = 2 X lU K , m - 

^00 Q W 

o 33 

2 X 10 cm , and M = 2 X 10 gm , we find that is 5.5 

o , • o 

approximately. Equations (13) and (14) then predict the pressure far 

-5 2 

from the sun should have the value 5 X 10 dyne/cra and the number 
7 -3 

density 2 X 10 cm . Since in the solution we have matched tem- 
perature in the interstellar medium, the density given by equation (14) 
being eight orders of magnitude higher than that measured in the inter- 
stellar medium implies that the pressures are out of balance by a like 
amount. (The pressure of the interstellar magnetic field being at best 
of the order of 10 dyne/cm offers no help in supporting the static 
corona.) The model thus grossly failing the condition for static equi- 
librium, we conclude that the only available mechanism for dropping the 
inexorable pressure of the solar corona is via action against the inertia 
forces of dynamic expansion. 
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Before considering in the next model the consequences of such 

expansion we close this discussion by noting that at* earth distance from 

5 o 

the sun the temperature given by equation (8) T = 4.34 X 10 K is 

€ 

comparable with though larger than the electron temperature 1,5 X 10 K 
measured by earth satellites (Montgomery, 1972), Since heat conduction 
is a process that is defined in the local rest frame of a gas, the 
process is not inhibited by dynamic expansion but in fact enhanced by 
the larger temperature gradients resulting from expansion. We there- 
fore expect heat conduction to play a significant role in any compre- 
hensive model of the solar corona. 


b. Dynamic , polytrope corona 

The model assumes steady, spherically symmetric flow with mass 
conservation leading to the first integral of the continuity equation 


the constant mass flux 
2 

qj -r p V r , 


/I 


where v is the bulk velocity and where p is the mass density. 


The equation of motion is 




4a _ 

dr 


The polytrope law 


(l.lb-2) 




where a is a free parameter is taken to replace the energy equation. 
Use of the law with equations (1) and (2) serves both to make a tract- 
able mathematical problem and a useful tool in analyzing the nature of 
solutions of the physical problem being modeled. 
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After dividing by the mass density p an»d substituting equation (3) , 
we integrate equation (2) over the inteirral r^ to r- yielding 


2 2 

I- Jk 

2 "2 


a p 


<y-L Of 


( 


c^-1 

P 



-G (i - i) . (l.lb-4) 

o 


We define at the base of the corona the characteristic speed 
2 

bv the relation c = ap /o and the nondime ns ional parameter 
^ o o o 

e ss which is seen to reduoe to the same expression as 

o p r 

o 

for the static corona. Thence making the changes of variable 


X = 


r 

o o 



permits recasting equation (4) in the caoonical form 



(l.lb-5) 


(l.lb~6) 


Equation (1) may be represented equivalently as 


y 



9 


2 

I c 
o o 


(l.lb-7) 


Since in equations (6) and (7) x is an explicit function of r, the 
problem Is to solve the equations for y snd z as explicit functions 

of X , 


The solution is obtained with the a£d of the substitution 



(1.1b -8) 
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Solving equations (7) and (8) for z we get 


z = 


2 2 

X 1 ^ "Si 

2 


M 


(l.lb-9) 


while combining equations (7) and (9) gives 


y = 


(y-l 2 

afl - o&l 


(1. lb-10) 


iJjhen equations (9) and (10) are substituted in equation (6) there results 
with some arranging of terms 


4 -2(<y-l) 

M . M 

2 Qf-i 


2 21^ 

/i.i +1° + 

^ x X ^ 2 ^ a-l/\ X / 


(1. lb-11) 


Equation (11) is of the form 


- 2(a-L) 

f (M) = \ g (X) . 


(1. lb-12) 


In equation (3) a = 1 describes an isothermal corona which does 
not admit finite heat conduction and therefore requires a heat source 
which is unidentifiable except near the sun. Since ot - y - 5/3 describes 
the adiabatic corona, i.e. no heating at all, we are led to investigate the 
topology of continuous solutions M(x) of equation (11) for a between 
1 and 5/3 . To clear up a point of confusion in Bondi *s related dis- 
cussion, in the solutions to be obtained for general a the characteristic 
1/2 

speed c = (np/p) is not the speed of sound a and accordingly, equa- 

tion (8) , M is not the Mach number M but these objects are related by 

c = a and M = . (1. lb-13) 
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To begin the analysis, since M is positive the left hand side of 
equation (11) f(M) is positive definite. The constant X (related to 
the mass flux) being positive, solutions of equation (11) require that 
the function g(x) be positive definite also. Since ~ approaches zero 
as r approaches infinity, solutions only exist when 



(1.1b -14) 


As we shall confirm by the solution the flow starts out very slowly in 

2 1 

relatively dense gas and y^/2 is a very small number compared to — 

Consequently with x = oi/e , from relation (5) , the condition (14) 

o o 

effectively requires 


Of < 



(1. lb-15) 


Since as we have seen in the discussion of the static corona model 

c has the value 5.5 approximately, relation (15) indicates that or 
o 

should be less than 1.22+ . 

Analysis of f (M) reveals that at M - 1 the function has a 
minimum of value 


f 

m 


0+1 

2(a~l) 


The function g(x) 


also has a minimum of value 



1 0+1 
^ 4 of-1 


4(1 + 


0-1 2 
— 5^0 


"o> 


5-3g 

Q+l 


5-3o 


0-1 

o 


(1. lb-17) 
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The minimum of the function g is located at 


5 - 30 ' 


X = 

m 




(1. lb-18) 


Of o 


Taking differentials of both sides of equation (12) provides the 
relation 


f dM - X 


cH-1 


g'dx . 


(1. lb-19) 


Equation (12) together with the differential relation (19) suggest two 
kinds of physically meaningful solutions to equation (11), 

The first type is the ‘’critical solution” obtained when \ is so 
chosen that the minima of f and g jointly satisfy equation (12). 

In this case f** and g^ both change from negative to positive through 
zero at the same point in the flow and the fluid accelerates smoothly 
through the critical point leading eventually to supersonic flow. The 
critical value of X is denoted by X^ and is given by 




C&l 

2(or-l) 



Q+1 

2(a-l) 


X 

ID 


-(3-3g) 

2(a-l) 


(1,1b- 20) 


Combining equations (20), (10), and (5) we have that the velocity at 
the critical point is 


V 

c 


2(1 yl . ^ 


1/2 


5-30f 


(1. lb-21) 
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Similarly with equation (9) the critical point density is 


p = p 


. 0^“^ 2 Q^-l . 

^ 2 ■' o a o ' 

5-3or 


O'-! 


(1. lb-22) 


The corresponding pressure follows from equation (2) as 


P ” P 
*^c o 


,,, a-1 2 a-l 

2(1 +— 7 ^ y T“ € ') 

^ 2 o a o' 

5-3a 


Of 

a-1 


(1. lb-23) 


Thence the perfect gas law gives for the temperature 


T = T 
c o 


c^-1 - a~l 

— V 

5-3or 


(iab-24) 


To exhibit the solution far from the sun, equation (11) approaches 
asymptotically 


_ 4 _ - 2 ( 0 - 1 ) 


M "*•’• = X 


/ 2 ^ 2 ‘'•o \ 

( - — j * 


2e_. ^ 


<rM 


(1. lb-25) 


Combining equations (5), (10), and (25) gives the asymptotic velocity 


V = c 
00 o 


2^2 

y + — ^ 
o or-1 


2e 


Of 


1/2 


(1. lb-26) 


which is seen to be a constant. 

The critical mass flux from equations (7) , (18) , and (20) is 

5-3Cg 
2(a-i)^2 

-iVo"! ‘ (l.lb-27) 


The density at large distance from the sun then varies as 


Offl 


% 


= (. 1 ) 


1 ) 


4(1+— ;5T- y €J 

2 •'o g Of 

5-3ce 
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'Pc -2 


p = r r 


(l.lb-28) 


Accordingly the pressure approaches asymptotically 


P = P, 




(1. lb-29) 


and the temperature varies like 


/tPcNa-l -2(a-l) 
T = T I-;; I r 

° 1^0 V 


(l.lb-30) 


The polytrope solution of the second kind, the subcritical solution, 
obtains when \ < \^ and the minimum in g is reached before f reaches 
minimum. In this case going through the derivative g' changes 

dM 

from negative to positive while f' remains negative. Consequently ^ 
changes from positive to negative and M is caused to retreat back along 
the lower branch of f from the maximum value M obtained as the solu- 
tion of 


— * 

f (M ) - 


OH-1 


-5+3g 

ot+1 


4(a-l) 


(1. lb-31) 


On the lower branch of f(M), M approaches zero as f(M) approaches 
Infinity, corresponding in this solution to g approaching infinity 
with X on the upper branch of g . Thus in the subcritical solution 
far from the sun equation (11) approaches asymptotically j 


- 2(g-l) -2(g-l) 


- 0+1 
M 


4(g-l) 

/ 1 + ail y 2 . Oil £ ) xa+i- 

^ ^ ^ 2 ^o 0^0^ 


(1. lb-32) 
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2 

Solving equation (32) for M and combining the results with equations 
(10) and (5) gives for the far field velocity. 





g~l 

a 




(1. lb-33) 


Inspection of equation (33) in the light of equations (7) and (5) shows 
that the density far from the sun approaches the constant limit 


Pco 




2 ■'o 


1 

q-1 > a-1 

"o> 


(1.1b -34) 


Accordingly from equation (3), the pressure is seen to approach the 
constant value 


P 


09 


Po + 


g-1 

2 


2 

3^0 


a 

g-l X g-1 

€ ) 

g o' 


(1. lb-35) 


Finally, the temperature also appracches a constant 


In summary of the character of the two kinds of solutions, the 
critical solution admits a flow that accelerates continuously approaching 
a constant supersonic velocity far from the sun. The density correspond- 
ingly decreases as r and the pressure and temperature both similarly 
approach zero far from the sun. Given conditions in the corona, the criti- 
cal solution for any given a uniquely defines a mass flux from the sun. 
The subcritical solution admits a flow that accelerates continuously 
to some maximum velocity less, than the speed of sound and thence slows con- 
tinuously approaching zero velocity far from the sun. The velocity history 
is carried out in such a manner that the density, pressure, and temperature 
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monotonically drop throughout, approaching constant miniroum values far from 
the sun. Choosing a value of a does not uniquely define the mass flux 
since a solution exists for all values of the mass flux less than that 
prescribed by the critical solution. Since the velocity far from the 
sun depends linearly on the mass flux and varies inversely with the square 
of the distance from the sun, the asymptotic values for the density, pres- 
sure, and temperature do not depend on the mass flux in the subcritical 
solution. 

Table 2 sets forth values for both types of solution for a = 1*15 
and a = 1,22 at salient points in the flows. The value a - 1.15 was 
chosen because it best reproduces in the critical solution the super- 
sonic velocity measured in the solar wind at earth distance from the 
sun. The value a - 1.22 was chosen because it best matches in the sub- 
critical solution the pressure in the local interstellar medium, 

Cwiivpariiig both the values of mass flux and velouiLy In cue ciiLiual 
solution at a = 1.15 and the subcritical solution at a = 1.22 


Table 2. Properties of tuo poly trope solutions for dynamic 
expansion of the solar corona. 


critical 


subcritical 


\ a 


r /R 
c o 


n 


1.15 

1.22 

• 


6.6 

2.15 



1.4 X 10"^ 

2.4 X 

10"^ 


, ■ ^12 -1 ^ -1 
1 X 10 gra s ster 

1.7 X 

o 

00 

-1 

s 

1.2 X 10^ cm s"^ 

2,1 X 

10^ cm 

-1 

s 


7 -1 

4 X 10 cm s 




0.03 cm 


-3 


-13 -2 

10 dyne cm 


0 



evidences the dramatic nature of the labels coined by Parker, "The 
Solar Wind”i and Chamberlain, "The Solar Breeze", in their now famous 
competitive advocacy of the respective solutions to describe the state 
of the solar corona (Parker, 1960 , 1963 ; Chamberlain, 1960 , 1961 ), 

In the perspective of conditions as measured by instrumented sa- 
tellites in the interplanetary region, Parker had the better of it, 
for from the coronal base to earth the solar wind undoubtedly follows 
a solution close to the critical polytrope. However,, at large distances 
from the sun the critical solution which predicts vanishing density, 
pressure, and temperature suffers a catastrophic failure and there a 
chamberlain- like solution is appropriate. 

Parker resolved this apparent paradox in a third ad hoc model which 
we take up in the next section. From the point of view of results de- 
veloped in the present section, the model provides a synthesis featuring 
discontinuous branching via a shock transrtron from a supersonic (critical) 
solution to a subsonic (subcritical) one. An additional feature of the 
third model providing greater physical content is that the sun is permitted 
to move relative to the interstellar gas. 
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2. Shock Termination of the Solar Wind 

The model to be presented follows the development brought forth by 
Parker (1961,1963) to demonstrate characteristic phenomena of the combined 
flow field resulting from the clash of the supersonic solar wind with the 
interstellar medium. As the result of certain bold idealizations of the 
physics of the problem made in the interest of obtaining a tractable 
mathematical formulation, Parker emphasized the "qualitative” nature of 
his results, disclaiming "quantitative and definitive" accuracy. 

The supersonic character of the flow admitting such a formulation, 
here and throughout the remainder of the thesis the discussion is from 
the point of view of an initial value problem beginning at 1 a,u. 

There the properties of the flow are now thoroughly established by ob- 
servation, (See for example reviews by Wolf, 1972; Gosling, 1972; 
Montgomery, 1972; and Schatten, 1972,) Typical numbers used in sub- 
sequent computations appear in Table 3, 


Table 3. Properties of the 

solar wind observed at 1 a.u. 

TEMPERATURE 

ELECTRONS 

1.5 

X 10^ 

PROTONS 

5 

X 10^ \ 

PARTICLE DENSITY 


-3 

ELECTRONS or PROTONS 

5 

cm 

MAGNETIC FIELD 
STRENGTH 

5 

-5 

X 10 Gauss 


Parker’s analysts is reviewed below with considerable amplification 
and one important extent ion. Here the material will serve the original 
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expository intent. But in later sections it will serve as an analytical 
tool in conjunction with theory to guide the development of and to ex- 
plore features of the solution to be expected from a considerably more 
detailed physical model whose construction and analysis are among the 
major subjects and purposes of this thesis. As the result of the inter- 
play useful light will be shed on the question of the accuracy available 
in the approximate solution of Parker's model as extended, the analytical 
solution being presently the only general description of the flowfield 
extant. Finally, wtiile not a task taken up in the thesis, the results 
of the analysis are expected to be useful in guiding the future develop- 
ment and application of numerical procedures to obtain a full solution for 
the flowfield described by the more exact model. 

Following ideas suggested in papers by Clauser (1960) and Weymann 
(1960) Parker constructed a solution for the termination of the solar 
wind by explicit ely assuming the existence of a shock transition. The 
shock provides a stable mechanism by which the very low pressure in the 
supersonic regime of the far interplanetary solar wind can be raised to 
match the much higher pressures of the local interstellar medium in the 
vicinity of the boundary where the two media must ultimately interact. 
Steady flow modeling average conditions is assumed in a non-rotating 
heliocentric reference frame, 

Parker made the following assumptions. The weakly magnetized gas 
of the interstellar medium is taken to be incompressible, implying the 
flow is subsonic in its passage of the sun. Also, the flow of the inter- 
stellar gas, termed "the interstellar wind", is assumed to be non-viscous, 
irrotational, and to form a sharp, mutually im.penetrable boundary with 
the magnetized post-shock solar wind. 
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To decouple the transonic flow problem governing the location, the 
geometry, and the conditions over the shock surface ‘from the solution 
of the ensuing subsonic flow problem in the solar wind, a spherical 
shock surface is assumed across which the uniform conditions on either 
side should obey the Rankine-Hugoniot normal shock relations. Since 
the Mach number of_the preshock supersonic solar wind is very high, 
the strong shock relations are invoked. Across a strong shock a highly 
supersonic flow becomes a moderately subsonic flow and the post shock 
solar wind can be treated as an incompressible fluid to a fair approxi- 
mation. Finally, the subsonic solar wind is assumed non-viscous and 
irrotational, also. 

With the above assumptions, both the post shock solar wind and the 
interstellar wind satisfy the relations of ideal fluid flow 


C X If ^V/ — > K — 


f't •> 1 \ 


curl (v) = 0 , and 


( 1 . 2 - 2 ) 


P Y~ *** P ~ constant . 


(l«2-3) 


Since the density p in either regime is constant the velocity 

1/2 _ 

V in equations (1) and (2) may be repleced by (p v = v") . The 
substitution results in the equivalent set of equations 


div (v") = 0 , 


curl (v') ~ 0 , and 





v". v' 


~ + p = constant 


( 1 . 2 - 6 ) 
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In these equations , no density explicitly appears and the constant of 
the Bernoulli equation is the same over all stream lines of both flow 
regimes, the constant being the stagnation pressure of the interstellar 
wind. Consequently there is nothing in the equations to distinguish 
them between e it tier flow regime and they represent the flow in both. 

The problem mathematically then has been reduced to that of a single 
ideal fluid of uv\lt density flowing in both regimes subject to certain 
boundary conditions at » and at the shock surface where the solar wind 
part of the flow field first becomes subsonic. 

To recover the solutions for the real fluids in either regime 
from the solution of the ideal fluid equations one simply applies the 
inverse transformation 



with appropriate density in corresponding regions of the ideal fluid 
solution. Since the pressure p in the Bernoulli equation (6) was not 
affected by the t: Vans format ion of the original equations, the pressures 

obtained in the Solution of the ideal fluid problem carry back unchanged 

to the respective regimes of the real fluid problem. 

Equations (A) and (5) admit a velocity potential that is the solution 
of Laplace’s equation. The solution of the equation is effected by 
superposition. 

At large distances from the sun (infinity) the effects of the solar 
wind are no longeT felt and the fluid appears as an incompressible uni- 
form parallel flow. At the spherical shock surface by assumption the 

solar wind flows Vadlally with uniform subsonic velocity and hence is 
indistinguishable from the flow of an incompressible point source 
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matching conditions at the distance of the shock radius. The flowfield 
of the combined subsonic flows is, then, to be constructed as the super- 
position of the appropriate parallel and source flows to which the com- 
bined flow is assumed asymptotic at the boundaries. By the vanishing 
nature of the source flow at infinity, there the solution so obtained 
will be exact. Since the parallel flow does not vanish identically at 
the shock surface, but depending on parameters may be small there com- 
pared to the source flow, the solution in the vicinity of the shock 
will only be an approximation more or less in keeping with the assumed 
sphericity of the shock surface and other idealizations of the problem. 

To Solve the subsonic flow problem as indicated above requires the 
location of nnd determination of conditions at the shock surface. This 
is accomplished using the integral of the continuity equation, with con- 
stant velocity in the supersonic region, to project conditions at the 
earth to the shock; using the strong normal shock relations to carry 
conditions across the shock; and using the incompressible Bernoulli 
equation to project conditions imposed by the interstellar medium back 
onto the shock to achieve a match up as detailed below. 

The Integral o£ the continuity equation in the preshock interplane 
tary solar wind gives 


2 ^2 2 
p V r = Pi R = Pg 



( 1 . 2 - 8 ) 


where subscript refers to conditions immediately on the supersonic 

side of the shock and subscript e refers to conditions in the solar . 
wind at earth distance from the sun. The radial distance of the shock 
surface from the sun is R . With the constancy of the velocity beyond 
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the orbit of earth, the density in the supersonic solar wind is seen 
to vary with the reciprocal of the square of the radial distance from 
the sun. 

The Rankine-Uugoniot normal shock relations are 


Pi - P2 


(1.2-9) 


Pi + Pi = P2 ^2 ^2 


( 1 . 2 - 10 ) 


2 y-l p, 2 y-1 p, 


( 1 . 2 - 11 ) 


where the subscript 2 refers to conditions immediately on the sub- 
sonic side of the shock. In the case of strong shocks » 1 , 

2 

?2 21 ^1 '^i Pi derived relations 


p2 7+1 
f>l ” 7-1 


apply 


( 1 . 2 - 12 ) 


Since for a monoatomic gas such as the fully dissociated and ionized 
hydrogen of the solar wind y - 5/3 , the right hand side of relation 

(12) has the value 4 . Substitution of relations (12) in (10) and 
dropping the negligible pressure gives for the post shock pressure 

in the solar wind 

P2= ^ Pi 
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Beyond the shock the incompressible Bernoulli equation holds in 
the subsonic solar wind so that immediately after the shock 

2 

^2 

^2 ^2 T " constant = (1.2-14) 

where is the pressure at the stagnation point. 

But the stagnation point pressure also obeys the incompressible 
Bernoulli equation in the interstellar wind region where 


a. 1 2 

= Pi + 2 '’l ''i 


1 2 
Pi- +2^1% 


(1.2-15) 


In equation (15) is the velocity of the sun with respect to the 

undisturbed interstellar medium and the pressure of the undisturbed 
medium is of the order 




8n 


+ p. a T 
1 


i 


(1.2-16) 


taking account both of the scalar pressure due to the interstellar 
magnetic field B^ and the thermodynamic pressure of the interstellar 
gas treated as atomic hydrogen. 

In the light of relations (15) and (16), the right hand side of 
relation (14) is a constant dependent only on the properties of the 
undisturbed interstellar medium and the motion of the sun through the 
medium. 

Substitution of relations (12) and (13) in the left hand side of 
(14) then gives 
2 

^5TPi T = Ps • 


(1.2-17) 
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Relation (17) is viewed as selecting the place in the supersonic solar 
wind flow where the momentum flux is such that if the shock occurs there, 
the post shock conditions will match up with a subsonic solar wind solu- 
tion the constant of whose Bernoulli equation is p as given by relation 

s 

(15) with (16). Finally, substituting for in relation (17) the 

value in terms of the parameters at earth and the radial distance from 
the sun through relation (8) we ge^ that the radial distance of the 
shock from the sun may be found as 


R = r 


(St) 


2 

P V 
e e 


2 P. + P.V 

lea 1 O 


1/2 


(1.2- 18) 


With R as given by relation (18), the immediate post shock 
velocity and density obtained from relations (12) and (8) are 



V 

1 Y^-V 

e 

= (^) 

2 

r 

e 

\Y~1/ 

2 


and 


(1-2-19) 


( 1 . 2 - 20 ) 


\ 


For the source flow representation of the post shock solar wind the 


source strength is given by 


Q = Att 


( 1 . 2 - 21 ) 


where v''' is the velocity of the source component of the ideal 
SR, 

fluid at the point R and is given by 


1/2 


sR 


= P. 


( 1 . 2 - 22 ) 
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Combining relations (19), (20), and (22) in (21) yields the source 
strength as a function of R . 


1/2 1/2 

4tt r R I P V 

e \Y+1/ e ( 


(1.2-23) 


J 

The potential function and Stokes stream function for the source com- 
ponent of the ideal fluid flow are respectively 


9. = “ 


_ 2 _ 

4ttt 


and 


I (1 - cos e) 


(1.2-24) 


(1.2-25 


where ^ is the flux in a cone of given semivertex angle 0 
s 

measured with the axis of flow in the direction of motion of the sun 
through the interstellar medium, i.e. in the direction upstream of 
the sun with respect to the interstellar wind. 

For the parallel flow representation of the far interstellar wind 
the velocity of the parallel flow component of the fictitious fluid is 



(1.2-26) 


The corresponding components of the potential function and Stokes 
stream function for the ideal fluid are 


cpp = - U X and (1.2-27) 

if = - TT U y^ 


(1.2-28) 
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In the equations, x is a cartesian vector component measured along 
the axis of flow, positive in the upstream direction, with respect to 
the interstellar wind, and y is orthogonal to x . 

Adding components for the source and parallel flows gives the 
potential function ' ^ 

tp = - U X - . (1.2^29) 

Similarly the Stokes stream function for the ideal fluid representation 
of the combined flows of the solar and interstellar winds is 

^ = 2 (1 » cos 0) - TTU y^ • (1.2-30) 


Differentiating equation (29) with "respect to x and setting the 
result to zero on the upstream axis of the flow places the stagnation 
point at 



(1.2r31) 


From relations (23) and (26) we find for the stagnation point location 


•k 

X 

o 



1/4 


P V 
e e 


Pi % 


(1.2-32) 


In the equation the asterisks signify quantities measured in astro- 
nomical units, i.e. normalized by r^ (See equation (18).) , 
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Combining equations (30) and (31) we find the streamlines in the 
X , y plane corresponding to the intersections with* streamsurfaces of 
constant velocity flux f are given by 

2 2 

y = 2x^ (1 - cos 0) 4* constant. 

The regions of the flowfield corresponding to the solar and inter- 
stellar winds are separated by the stagnation point streamline obtained 
from equation (33) with the constant set equal to zero. Streamlines in 
the interstellar wind regime correspond to positive constants and those 
in the solar wind regime to negative constants. The most negative con- 
stant yielding a real streamline is (-4 x^) corresponding to the axial 
streamline in the downwind direction. At large distances downstream of 
the sun the streamlines become asymptotically parallel corresponding to 
circular cylindrical stream surfaces. In particular tne radius of the 
stagnation point stream surface approaches asymptotically 

y == 2x , (1.2-34) 

*^o o 

With equation (31), equation (29) may be represented as 
2 

cp - . U (X ) . (1.2-35) 

Then the ideal fluid velocity components obtained by differentiating 
equation (35) are 


(1.2-33) 
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/ 2 

/X X 

V-: = u - 1 


v; = u 


and 


(1.2-36 


(1.2-37 


Finally, making use of equation (7) and^ (26) we have the real 
fluid velocity components 


V = V 


\ 1/2/ 2 

Pi\ ho ^ 


o \ o 


V = V 


V 1/2 2 


o\ p 


- 1 ) and 


(1.2-38) 


(1.2-39) 


where setting p = p. gives the velocity in the interstellar medium and 

/r \2 ^ 


P == P- 


gives the velocity in the subsonic solar wind. 


2 y-l / ■ e 

Equation (38) predicts that at large distances downstream of the sun 
the velocity in the subsonic solar wind should approach the constant 
value 


= V , , , 

X o Xy+l p 


\y+l pJ 


1/2 


R 


(1.2-40) 


The speed with which the sun moves relative to the local inter- 
stellar gas is presently quite uncertain. Astronomical measurements 
put the speed of the sun with respect to the nearby stars at 20 km/s 
(Allen, 1963), In the absence of other better indications, this value 
has been the standard interstellar wind speed for purposes of cajlcu- 
lations and physical inquiry. 

However, consistent with both the expanded knowledge of the temper- 
ature and density in the interstellar medium and the' developing model 
for the interaction of the solar wind with interstellar neutral hydro- 
gen, recent interpretation of measurements of interplanetary Lyman - cx 



radiation provide both grounds (Bertaux and Blamont, 1972) and a 
method (Thomas 1972) for making a new, more direct estimate. In 

the reference cited, Thomas finds the speed of the interstellar wind 
to be on the order of 6 km/s . 

Since the model used by Thomas is somewhat incomplete, for the 
present, we regard the value given as indicative. The subject of 
Lyman - cx scattering with inferences to be drawn from it is discussed 
further in section 4, and again in Chapter m, in conjunction with cal- 
culations of interstellar wind speed based on our own model. 

We now explore the possibility of using the foregoing solution of 
Parker to provide a rough approximate description of the combined 
solar wind - interstellar wind flow field under conditions typifying 
current estimates. To point up a separate major difficulty in the 
model as originally formulated and directly interpreted we initially 
push the asSuiViption of incompiessIbiliLy in the iuLersteilar mediura 
to the limit by assuming the sun to move with sonic speed ( ~ 16 km/s ) 
relative to the ionized component of the interstellar gas. Later, 
with the model reinterpreted when we reduce the assumed interstellar 
wind speeds toward the substantially subsonic estimate of Thomas, rea- 
sonable consistency is achieved in the incompressibility assumption. 

To proceed accordingly, we choose from Table I the particle density 

-3 

of the interstellar neutral gas to be 0.1 cm , the density of the ionized 

-3 

fraction to be 0.01 cm , and the interstellar magnetic field strength 
“6 

to be 3 X 10 gauss. Putting numbers in equation (18) where and p^^ 

include contributions of the ionized as well as neutral interstellar gas, 
we find 125 a.u, as the estimate of the distance to the shock termination. 

If the same procedure is followed for equation (32) using the foregoing result 
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for R , the stagnation point is found to be at 100 a.u., well within 
the shock radius. The- latter is an obviously absurd- result from which 
we infer that some assumption of Parker *s model must be grossly incon- 
sistent with the physics of the problem to which the model is being 
applied. Understanding of the difficulty is promoted through a relation 

for X /R obtained from equations (18) and (32) . ^ 

o 

The relation is 


R 


2(y-l ) P j.q* t 1/2 


y+3 


p. V 

1 o 


2 T 


1/4 




(1.2-41) 


Viewed as representing an ordinary gas dynamics problem, equation (40) 
may be recast 


'izL li i, + 1] 1 
7+3 \y / J 


1/4 


(1.2-42) 


With 7 given the value 5/3 , the coefficients and 2/7 have 
the values 1/7 and ‘6 /5 respectively. Hence in order^that the model 
place the stagnation point outside the shock surface > 1) , 


equation (41) requires that < -^ or the free stream Mach number 

M < 0.45 c 

CO ^ 

Based on Thomas’ estimate of v , the Mach number of the sun 

o 

through the neutral interstellar gas is 0.5 . Quite evidently then 
Parker’s model in the purely gas dynamic interpretation of the problem 
predicts that in the upstream direction the assumed incompressible 
source flow of the solar wind will be deflected back and around the 
flow of the interstellar gas before the source flow ever reaches the 
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place where the shock in the real supersonic solar wind ought to be. 
Parker himself noted (1961) that the representation should apply 
only for very slow passage of the sun through the interstellar gas. 
Somewhat ironically but fortunately for the vitality of his 
model, in fact, the bulk of the interstellar wind does not flow as 
Parker originally hypothesized. As recognized by Holzer and Axford 
(e.g., Axford, 1972) and as we shall see in more detail in the next 
section where we discuss particle interactions, outside the shock the 
mean free path for a neutral hydrogen atom to interact with anything 
is much greater than the characteristic dimension of the flowfield, as 
measured by R for example. Consequently the neutral hydrogen is 
not deflected as a continuum fluid forming a definable boundary with 
and flowing around the solar wind in accordance with Parker *s model. 

But the numerically much weaker ionized component of the interstellar 
wind being a magnetized plasma of short gyro-radius interacting with 
another similar plasma does form such a boundary and is deflected to 
flow in such a way. To the zeroth order approximation of the present 
discussion, then, the neutral hydrogen is regarded as passing without 
effect through both the ionized interstellar wind and the subsonic 
solar wind. 

If we now examine the assumptions of Parker *s model in the ‘light 

* 

of the foregoing, we note that the relations for the location of the 
shock distance are a statement about exchanges of momentum between 
the two media. As the result of momentum conservation in the x direc- 
tion in single particle collisions, it does not appear that on the axis 
in the forward direction it should make too much difference whether a 
fraction of the overall momentum flux of the interstellar wind is 
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exchanged with the solar wind within the shock as in the case of the neutral 
hydrogen or whether it is all stopped at a stagnation point as in Parker's 
assumption. 

On the other hand, given the location of the shock the location of 
the stagnation point in Parker's model can be seen to derive from kine- 
matical relations and is a statement about the conservation of velocity 
flux in the flow of two incompressible, irrotational fluids that must 
evade each other. Evidently the two such fluids here should be the 
solar wind, in principle augmented by the ionization of neutral hydro- 
gen transported into the solar wind, and the ionized interstellar wind. 

Consequently while in equation (18) it is reasonable to take 
as the sum of the mass densities of the neutral and ionized components 
of the interstellar wind, in equation (32) should reflect only the 

density of the ionized component. When this is done in equation (32) 
the prediction is that the solar wind be bounded at 185 a.u,, a dis- 
Lciiictr well ouLblde the siiock. 

Following the corrected procedure but with Thomas* estimate of 
the interstellar wind speed (6 kra/s) we find 145 a.u. for the shock 
radius and 328 a.u, for the distance to the stagnation point. Stream- 
lines of the flow field are sketched in Figure 1. 

An indication of the internal consistency of the model as extended 
is available from equation (38) evaluated on the upstream axis of flow 
in the subsonic solar wind. With equation (32) and x and r taken 
equal to R , it may be seen that the first term of equation (38) 
gives the post shock velocity, equation (19); and accordingly the second 
term represents the error at the shock. The relative error in the sub- 
sonic solar wind velocity in the vicinity of the shock is hence measured 

by R^/x^ , about 20% for v =6km/s , 

o o 
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Using the value for the shock radius 145 a.u*, we see from equation 
(40) that the velocity in the subsonic solar wind downstream of the sun 
is of the order of 3 v^ or 2 x 10 cm/sec. We recall the speed 
of the sun to be equivalent to about Mach 0.5 in the interstellar gas. 
Thus downstream of the sun the flow of the solar wind has the character 
of. a transonic jet in the interstellar medium. ^ 
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3, Heating and Flowfield Temperature 

We now complete the gross characterization of the flowfield beyond 
earth orbit with a discussion of the temperature. The extended shock 
termination solution provides the general description of the flowfield 
on which the analysis is based. 

In the supersonic regime of the solar wind flow the temperature is 
governed by a competition between heating by conduction and cooling by 
expansion. Use can be made of the heat conduction and critical poly- 
trope solutions examined in the previous section to get a qualitative 

view of the situation. The polytrope solution predicts constant velocity 

-2 

and resulting density and temperature functions that vary as r and 

-2(a-l) , 

r ' ' respectively. 

The volume rate of heating Q associated with any prescribed flow 
can be got from a representation of the first law of thermodynamics, 
rhe hear equation for steady flow 


Q = p div V + p V . grad e 


(1.3-1) 


Expressing div v in the radial coordinate with spherical symmetry, 

the internal energy e by (p/p) > P the polytrope law 

P ” P.»(p/p^) invoking the critical polytrope solution yields the 

12 C 

result 


Q « 


2p v 
*^e e 


(3)(^) 


2crt-l 


(1.3-2) 


Evidently the requirement for heating is substantially reduced both with 
increasing a and distance from the sun. For a “ Y > the adiabatic 
solution, there is no heating. 
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The heating available from conduction is given by 

A 1 d , 2 dT- /I Q ON 

Q =:- div q = -"2 ^ * (1.3-3) 

r 

Neglecting for the moment the effect of the interplanetary magnetic 

field on conductivity in using equation (1.1-4) and again assuming 

the polytrope dependence for T leads to the result 

7/2 7<y-5 

« 2AT r 

Q = 1 (cv-1) (7<y-8) (^) (1.3-4) 

r 

e 

Inspection of equation (4) reveals no heating (Q > 0) is avail- 
able from conduction when 1 < or < 8/7 • When the upper limit of 

the range is substituted in r the temperature dependence of 

-2/7 

the accorded polytrope solution is found to be the r variation 

of the constant heat tiux solution of Chapman. 

Solutions in the range 8/7 < a < y - 5/3 can derive heating 

from conduction. It is interesting to note that the effective polytrope 
index best representing the flow of the solar wind from the base of the 
corona to earth a = 1.15 lies rather closer to the lower end of the 
range indicating the system adopts gradual heating of the flow from 
slow attenuation of the conduction heat flux. Indeed if one equates 
the expressions for heat supply and demand, equations (2) and (4), to 
get the exponents of r in consonance one finds 1.2 , furthering 

the conjecture. The transport coefficient of thermal conductivity A 
required to bring the coefficients of the two expressions to equality 
is found to be 5,6 x 10 ^c.g. s, units, a value consistent with theoretical 


estimates 
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In the vicinity of the solar poles where the interplanetary mag- 

* • 

netlc field lines make only small angles with the radial direction, the 
state of the system should be essentially as we have just described. But 
near the solar equatorial plane the situation is more complicated. 

There it is now well established the magnetic field lines follow 
an Archimedes spiral approximately according to the law 


tan 


-1 ’^e 


(1-3-5) 


where Q is the angle the field direction makes with the azimuthal 
direction. 

As presently understood theoretically and as we shall see in the 

■next chapter, the consequences of the magnetic field becoming more 

2 

azimuthal is to rapidly lower (by .sin C) the effective coefficient of 

thermal conductivity in the radial direction. Beyond several earth - 

radii from the sun the heat flux is effectively blocked. Qualitatively 

- 4 

the. result is a temperature profile that falls off faster than r * 

(for ot - 1,2) which was found in the absence of the field. 

The effect of cutting off the heat flux can be partially analyzed 
through the dynamical form of the heat equation 


div 


2 GM 

P v(|- + ^ p/p - 


div q 


(1.3-6) 


Assuming spherical symmetry and conservation of mass (equation l,lb-l) 
we get from equation (6) the integral 


h- 

v-l 


R(T - Tp ^ (V* 


, GM r 

v ) + (1 - ) 

e' r ' TC ' 

e 




(1-3-7) 
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Equation (7) is seen to embody a prediction of the temperature in the 
supersonic region provided consistent estimates are available for the 
second and last terms. 

Since at earth the angle which the interplanetary magnetic field 

makes with the radial direction has only reached about 45 , heat 

conduction from the base of the corona has not yet been greatly effected. 

Hence for purposes of estimating the heat conduction flux at earth we 

-.4 

take the temperature profile locally to have the r dependence found 
in the absence of the field. Then taking account of the local inhibition 
of conductivity by the field the heat flux is given by 

T 

q = .4 sin^ ^ A , (1,3-8) 

e r 

e 

Taking T^ to be the electron temperature 1,5 x 10^ measured 

13 

at 1 a. u. gives the value 3.12 x 10 erg/gm for the right hand 

side of equation (7). Out near the shock (r 125 a. u,) the gravi- 
tational term has the value .88 x 10 erg/gm . Taking T in the 

5 o 

first term of the left hand side to have the average value 1, X 10 K 

of the electron and proton temperatures measured at 1 a. u. , one finds 

the term T^ = 4.16 x 10^^ erg/gm * 

2 14 

The term v /2 has the value 8 x 10 erg/gm which is substan- 
e 

tially larger than the others. Hence if we are to use equation (7) to 
estimate the temperature we must be able to gauge the dynamics rather 
closely. The extent to which the latter may be accomplished depends on 
the importance of heating to the remainder of the dynamical trajectory. 
To investigate the question we invoke the one dimensional compressible 
flow equation (Chapter II) in radial symmetry. 
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2 2 

(a ..y ) ^ 

V dr 


o 2 GM A 

% + -| + (v- 1 ) ^ 


(1.3-9) 


pv 


Using the relation for the speed of sound 


a « 2 7 a T 


(1.3-10) 


2 13 2 —2 ^ 

we find at earth = 2,77 x 10 cm s which is much less than v . 

® 2 

d V 

Consequently the left hand side of equation (9) is - (-^ ) and making 

use of the constant mass flux relation we write 


d 2a^ , ,, Q r 

d? <2 ) = - - T- - ^ ; 


(1.3-11) 


-1 -1 

At earth the first term on the r.h.s. has the value 3,7 erg gm cm 
compared to - .59 erg gm”^ cra**^ for the gravitational term. 

Beyond earth the thermal conductivity diminishing factor sin Q 
has the approximate dependence 


sin 


‘"ft) 


‘1-"^ c 1 


r ^ r. 


tan^ C I r > r 
’ e 


(1,3-12) 


Using the r functional dependence of equation (12) and again assuming 
the polytrope temperature dependence we find 

2 


Q = 


2(a-l)(7(y-6) A T 


7/2 


m 


la-3 


sin Q 
2 

tan Q 


(1.3-13) 


for the volume rate of heating by conduction. Comparing equations (4) 
and (13) we see that the relative effect of the interplanetary magnetic 
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field is to increase heating near earth (a ~ 1*2) but ‘rapidly to 

. J/ 

2 2 > 
diminish heating (by ), farther away. * . 

At earth the heating term on the right hand side of equation (11) 
then has the value -3,33 erg gm ^ cm ^ which is comparable with the 
value of the first term on the right hand side. Hence as a result of 
the closing magnetic field lines in the solar equatorial plane heating 
is important, in a relative sense, in the dynamics of the supersonic 
solar wind near earth. However substituting equations (13) and (10), 
with the polytrope law, in equation (11) we see that the terms on the 
right hand side vary with r respectively as r r ^ , r 

Then beyond earth, where ^ assumed, the heating term diminishes 

-2 

with respect to the first at least .as fast as r , Hence for purposes 
of estimating from equation (7) the temperature out near the shock 
(r ^ 145 a. u,) we can get a fair estimate of the integral of equation 
(11) by neglecting heating. 

The integral is formally available to us in the adiabatic poly trope 
solution (of =5 y = 5/3) already developed' (section 1,1b). The solu- 
tion is given relevance to the present problem by basing the initial 
point parameters ’ ^o * ^o conditions measured at 

earth. 


Since for the adiabatic solution c is the speed of sound a , 

e ^ e 

y = M = M the Mach number which is 7.7 at earth. The Mach number 
•^e e e 

being greater than one, the solution is on the ascending branch ^!of f(M) . 

The initial value for the non-dimensional distance x has the 

® / 

value 3^12 which is greater than zero, the location of the minimum 
of g(x) for the adiabatic solution. Hence is on the ascending 

branch of g(x) , also. 
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The initial values for x and M both being on ascending branches 
of their respective function curves, the solution we seek to equation 
(1, lb-11) belongs to the class of critical solutions. Since out near 
the shock x and M have values of the order of 400 and 20 respec- 
tively, the asymptotic formulas are appropriate. Then the result we 

2 2 2 

seek, equation (1, lb-26) gives v /2 = 8.33 x 10 cm s for v = 

7 -1 

4,08 X 10 cm s , an increase in velocity of only 2% beyond earth. 

Since by equation (11) the effect of heating a supersonic flow is to 
retard the velocity, the result given is an upper bound. 

13 

Thus the second term of equation (7) has the value 3,3 X 10 
erg gra ^cm ^ , Accumulating numbers found for the other terms pre- 
viously, we find 74,000 for the estimate of the temperature at 
the shock. But the polytrope law even with a = 1.14, i.e. unattenu- 
ated heat conduction, gives a temperature of only 25,000 ^K; and for 
a = 1,2, gives 14,000 ^K. Evidenciy the estimate of the temperature 
that we have derived from equation (7) is much too high. 

The source of difficulty lies not with the equation itself, 
which is fundamental, or with the estimate of the terras on the left 
hand side, the least well known one being an upper bound. Hence the 
problem must lie with the right hand side of the equation. 

Two possibilities suggest themselves. The first is that the 

-7 

classical transport coefficient of heat conduction 6 x 10 c.g.s. 

-2 -1 -1 

units (and the resulting estimate of heat flux 1,04 x 10 erg cm s 
at earth) is too large. The second is that the hypothesis of effec- 
tive extinction of the heat conduction flux at large distances beyond 
earth by the closing magnetic field lines is tod severe. 

That the former is true has been inferred previously on a variety 
of different empirical and solar wind- theoretical grounds by Montgomery 



(1972), Scudder (1972), Barnes and Hart le (1972), Whang (1972), and 

Cupennan and Harten (1972). Estimates of the required reduction based 

on solar wind theory range from a factor of 2 to 5.5 * Numbers at 

the high end of the range are consistent both with the predictions 

of the polytrope law and Scudder *s interpretation of the conduction 

heat flux in the steady solar wind derived from measured electron 

distributions. Corresponding in our calculation to a reduction factor 

-3 -2 -1 

of five, the heat flux 2 x 10 erg cm s is typical of the numbers 
quoted by Scudder in the ’’convergence region", i.e. the third moment of 
the distribution function cut off between 200 and 300 ev. (With re- 
gard to the cut off point, it is noted that the r.m.s. value of the 
Maxmillian distribution for the electron temperature at 1 a. u. is AO ev. 
If using the quoted reduced heat flux at earth, we recompute from equa- 
tion (7) the estimate of the temperature at the shock we find 14,000 ^K, 
the prediction of the polytrope law for a - 1.2 . 

In the reference cited, Scudder also has observations relevant to 
the possibility that the effective coefficient of heat conduction trans- 
verse to the interplanetary magnetic field is not greatly reduced with 
respect to the coefficient in the field free or field parallel case. 
Looking at the directional distribution of the heat flux vector about 
the magnetic field lines, he finds that the flux due to unsteady injec- 
tions of high temperature electrons into the solar wind is ducted by 
about a factor of three by the magnetic field lines. But for the steady 
electron population in the convergence region the direct interpretation 
of the data presented is that conduction transverse to the field is not 
greatly reduced, apparently only by a factor of about 0.7 , with re- 
spect to the field aligned component. 
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If in addition to the reduced heat flux at earth we accept the 
direct interpretation of Scudder’s data, namely that the steady com- 
ponent of the electron heat flux transverse to the magnetic field is 
transmitted with only a factor of 0«3 attenuation, the estimate of 
the temperature out near the shock is reduced to only 400 K « We 
take this value to be a lower bound on the temperature and the previous 
estimate 14,000 to be an upper bound. It is noted that the tem- 
perature predicted by the strictly adiabatic polytrope expansion is 
130 \ . 

Thus in the present state of uncertainty the range of temperature 
estimates in the supersonic solar wind near the shock essentially dupli- 
cates that 'in the interstellar medium. Hence we can take as a rule of ■ 
thumb that the temperature in the far supersonic solar wind is compar- 
able with that in the interstellar medium* Accordingly, the speed of 
sound is at most 20 km s^‘ and the pre shock Mach number is at least 
of the order of 20 . 

Because the Mach number in the preshock supersonic solar wind is so 
high, the post shock pressure and temperature are only very weakly de- 
pendent on the preshock pressure and temperature through the strong 
shock relations. Making use of the relations developed in section 1.2 
we* find that the post shock and stagnation pressures are given respec- 
tively by 
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The coefficients in equations (14) and (15) differ only by about 
Then the subsonic solar wind is only slightly compressible, the post 
shock density being given by equation (1,2-20). Hence the temperature 
in the subsonic solar wind may be considered isothermal with value given 

by 

2 

To == « 1-8 X 10*^^ \ . (1.3-16) 

2 (y^I) 

Thus we close this section by noting as a second rule of thumb that 
the temperature in the subsonic solar wind replicates that at the base 


of the corona. 
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4. Modern Developments - Particle Interactions 

In the historical context, one can equate the early development of 
the subject with the establishment of the flowfield by classical continuum 
methods, the focus of our discussion to this point. This phase spanned, 
six years culminating in the publication of Parker’s (1963) book. 

Signaled by the publication in the same year of papers by Axford, 
Desslcr, and Gottlieb (1963) and Patterson, Johnson, and Hanson (1963), 
the modern development is characterized by the occupation with filling 
in a wealth of physical detail missing in the early macroscopic models. 

Much of the effort has been concerned with identifying the important 
microscopic physical processes underlying bulk behavior within and between 
the solar wind and the interstellar medium. The progress of research, 
involving of practical necessity the rather piecemeal testing of concepts 
in the evolving corporate model, has now led to the point of understanding 
v/here it is timely to enter a new era of analysis from synthesis. 

Thus the present section summarizes, unifies and in some particulars 
clarifies the contributions to the physical model made in the recent liter- 
ature. The remainder of the thesis is then concerned with producing a 
coherent mathematical model adequate to explore features of the combined 
flow problem with the increased accuracy and for the subtlety of detail 
the added physics portends. 

The paper by Axford, a]L, (1963) touched on many noteworthy topics 
including the role of the Interplanetary magnetic field in the dynamics of 
the solar wind, a subject we treat in Chapters II, and IV, But the point of 
greatest significance to the present development, at least^ is the recogni- 
tion of the compelling likelihood of the charge exchange interaction between 
solar wind protons and atoms of neutral hydrogen. The neutral hydrogen 


is of interstellar origin and is able to invade the cavity occupied by 
the solar wind in the galactic magnetic field. The authors also ten- 
tatively estimated the effect of photoionization of atomic hydrogen by 
solar ultra-violet radiation in the Lyman continuum, 

Patterson, et al, then introduced the important notion of using 
measurements of extraterrestrial Lyman-a radiation to infer both pro- 
perties of the interstellar medium and structure of the solar wind through 
correspondence with predictions of an interdependent physical model. As 
the result of large wavelength dispersion, the Lyman-or was considered 
to originate in backscatter of solar radiation from hot interplanetary 
neutral hydrogen. The hot interplanetary hydrogen, in turn, was sup- 
posed to be the product of charge exchange reactions between hot. solar, 
wind protons and interstellar hydrogen which at that time was thought 
to be cold ( «100 °K) . Following Axford, et al,, the reactions were 
assumed to take place in a thin shell subsonic region. 

Hundhausen (1968) recognized the thin shell assumption to be falla- 
cious on grounds of the large mean free path for charge exchange. In a 
more systematic model, he then showed that the assumption the Lyman-of 
scattering should come from secondary hydrogen leads to an unreconcilable 
theoretical contradiction in the location of the shock surface in the 
solar wind. 

Semar (1970), following the work of Biermann, Brosowski, and Schmidt 
(1967) for cometary flow, pioneered numerical calculations of sjlar wind 
flow using one dimensional gas dynamic equations. The equations featured 
source terms for mass, momentum, and energy additions based on photo- 
ionization and charge exchange reactions with interstellar hydrogen. 



Blum and Fahr (1970,1971), Fahr 1972, Holzer and Axford (1970 a,b) 
Holzer (1972) and others have elaborated various aspects of the transport 
of interstellar neutral hydrogen into the supersonic solar wind region, 
the influence of charge exchange with this hydrogen on the location of 
the shock, and the interpretation of Lyraan-of scattering from such hydro- 
gen. (See Axford (1972) for additional references and a comprehensive 
review.) In the present context, all the forementioned work has suffered 
from the failing to incorporate both (1) a high temperature velocity disr- 
tribution and (2) the effect (see below) of solar Lyman-a radiation pres- 
sure on the interstellar hydrogen gas. 

Finally, Thomas (1972) brought to the field in modern context (but 
see also Brandt, .1960) the idea developed by Wilson (1960) that the radi- 
ation pressure of solar Lyman-a approximately cancels solar gravitation 
with respect to the transport of interstellar neutral hydrogen. Also in- 
corporating an isotropically distributed velocity component to roughly 
simulate a temperature distribution for the hydrogen, Thomas was able to 
gain an estimate of the speed of the interstellar wind. 

Using Thomas* estimate of the speed ( ^6 km/s) and the resulting 
estimates of the flowfield parameters based on the extended Parker solu- 
tion (sections 2. and 3.), we now turn our attention to the physical 
processes mentioned above and to demonstrating the scale of their effects. 
Based on results obtained in the earlier sections, Table 4 collects nomi- 
nal values for parameters of the solar wind regions and the interstellar 
wind . 

a. Photoionization 

A hydrogen atom exposed to solar ultraviolet radiation of wave- 
length at or below the Lyman edge 912 A may absorb a photon and become 
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ionized. The process thus contributes one free proton and one free 
electron to the plasma of the flow regime in which the event occurs. 

If the absorbed photon has a wavelength less than 912 A , correspon- 
ding to the hydrogen ionization threshold of 13.59 ev, the residual 
energy of the photon appears largely as kinetic energy of the freed 
electron. 

The photoionization rate per atom of hydrogen is given by 


or 


P 


J* f(v) tTp(v) dv 


(1.4-1) 


1 


Table 4. 

Nominal values for 
interstellar wind 

parameters of solar 
flow regimes 

wind and 

SOLAR WIND 


INTERSTELLAR WIND 

SUPERSONIC 

SUBSONIC 

IONIZED 

NEUTRAL 

NUMBER 

5 

0.001 

0.01 

0.1 

DENSITY^. 

r2 




(cm" ) 





TEMPERATURE 

10^ 

loV 

lo" 

10^ 

(°K) 





velocity 

4xlo’ 

< 10^ 

6x10^ 

6x10^ 

(cm/s) 





MACH 

> 7 

< 0.5 

0.4 

0.4 

NUMBER • 





— distances 

r in astronomical 

units. 



where v is the frequency at the Lyman edge, f(v) is the photon 
o . 

flux density, and a (v) is photoionization cross section. As the 
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result of the approximate v frequency dependence of the cross section 

above the Lyman edge, the integral is rapidly convergent; and the effecr 

tive wavelength band Is approximately to 2v^ • Using the solar 

u. V * flux measurements of Hintegregger , et al. (1965) ^ Banks and 

Kockarts (1971) have computed the photoionization rate to be 

-7-1 

1.5 X 10 s. at 1 a. u. The integrated flux over the effective 

wavelength band is 2.8 X 10^^ photons cm ^s ^ and the average photo- 

— -18 2 

ionization cross section a is 5.4 x 10 cm (Holzer, 1972) « 

-3 

Assuming the v dependence for the cross section and uniform 
photon flux permits an easy estimate of the average kinetic energy of 
the photoelectrons, i.e, the average excess energy of the absorbed pho- 
tons. When the integration is carried out over the interval v to 2v , 

o o 

the average photon energy is 4/3 hv^ where hv^ is the photon energy 
13.59 ev at the Lyman edge. 

Consequently the average energy residing in the photoelectrons is 

roughly 1/3 hv or 4.5 ev. This number is only half the 10 ev figure 
_o 

adopted by Semar from Biermann*s (1967) estimate for heavy molecules of 
cometary flow. Evidently the latter number is inappropraite for hydrogen 
The mean free path for photons is 


(N cTp)'^ (1.4-2) 

where N is the number density for hydrogen. At most, N has Ithe 

interstellar number density N = 0.1 cm and, consequently L 

o • / ^ 

18 5 / 

equals or exceeds 1.85 X 10 cm or roughly 10 a. u. Since the 

-3 

characteristic dimensions of the flowfield are less than 10 a.u,, 

attenuation of the ionizing flux is not a factor in the problem. With 
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the assumption of spherical symmetry, the flux then varies as r * 
Hence from equation (1)”, a niay be represented in tierms of the com- 
puted value (X at earth by 



(1.4-3) 


The mean free time t for hydrogen before photoionization is 

P 

given by the reciprocal of • The characteristic speed for the trans- ^ 

port of interstellar neutral hydrogen into the solar vind is the mean 
thermal speed at 10^ . Then the mean free path against photoioni- 

zation is the product of the characteristic speed 1,5 x 10 cm/s with 

•jf ^2 

the mean free time, yielding X - 0.6 r where as before the asterisks 

P 

signify quantities measured in astronomical units. 

The differential relationship between hydrogen number density and 
mean free patn is 


dN = - ^ N , (1.4-4) 

Assuming for purposes of the present rough estimation that the hydrogen 
moves radially toward the sun we have d^ = - dr . Then with the approxi- 
mate representation for the mean free path found above, equation (4) may 
be integrated to give N = exp (-1.5/r ) . Evidently a fraction on 
the order of 15% of the interstellar hydrogen is unable to survive photo- 
ionization to reach 10 a, u. while only 1%. suffers photoionization in 
passing through the subsonic region to reach the (shock) boundary of the 
supersonic region. 

The protons which are added thusly to the supersonic solar wind carry 
mass and, in conjunction with charge exchange, act through momentum 



conservation to retard the velocity of the far supersonic solar wind. 

A needed parameter in a subsequent calculation of velocity retardation 
(subsection c.) is the amount by which the proton flux of the super- 
sonic solar wind is augmented by photoionization. 

The conservation law form of the equation governing proton number 
density n in the steady solar wind (Chapter II.) is 

div (nv) = ;p 

where n is the volume rate of photoionizations given by 


n 

P 


o N = 
P 


2 

a r 

ep e 


N 

2 

r 


(1.4-6) 


With the assumption of spherically symmetric radial flow in the supersonic 
region, equations (5) and (6) may be expressed as 


I 

2 

r 


d 

dr 


(n V 




a 


ep 



(1.4-7) 


Recognizing the quantity under the differential operator to be the 
proton flux cp , we integrate equation (7) giving 


% 


a 


6£ 


e e 


r 


J N d r 
r 

e 


a N 
ep Q 

n V 
e e 


(r - r^) 


(1.4-8) 


b. Charge exchange 

In the charge exchange reaction, a proton is passage of a hydrogen 
atom acquires the temporally shared electron. Since in the case of charge 
exchange between protons and hydrogen the product species are identical 



in kind with the reactants, the process is said to be resonant. That is, 
the reaction requires no exchange of kinetic for cheraical potential energy. 
As a consequence the reaction is more probable the longer the reactants 
are in proximity; and, accordingly, the cross sections are higher for 
lower relative velocities between the reacting particles. 

The charge exchange cross section has been calculated over the energy 
range of our problem, 1 ev to 1 Kev by Dalgarno and Yadav (1953). 

At 1 ev, corresponding approximately to the case of protons and hydrogen 
interacting in the interstellar wind at 10 K , Dalgarno and Yadav find 
the cross section to be 47.3 X 10 cm . At 1 Kev, corresponding 
approximately to the case of supersonic solar wind protons colliding 
with interstellar neutral hydrogen, the calculated cross section is 
16,5 X 10 cm 

Over the subrange 20 ev to 1 Kev the charge exchange cross section 
has been determined experxmencaliy by kite, Smith and ijtebbings (196^). 

The measurements are found to be about 10% higher but qualitatively fit 
the calculated cross sections very well. Indeed to accuracy well within 
that of the experimental measurements, Fite, find the measure- 

ments fit the semi-empirical law for symmetric resonance cross sections 
(Dalgarno, 1957) « 

= 7.6 - 1,06 log^pE (l,4b-l) 

1/2 -8 . « . . 
where a ' has units of 10 cm and E is xn ev. 
c 

Expressing E in terras of the relative velocity of the particles, 
we find by squaring equation (1) the relation for the charge exchange 
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cross section 


a = 62.3 
c 


- 33.5 log^Q V + 4.5 



V 


(1.4b-2) 


“16 2 

In equation (2) the cross section has units of 10 cm and the velo 
city is in units of lO^cra/s. All charge exchange cross sections used 
in the thesis are derived from equation (2). 


The charge exchange cross section plays its role in the reaction 
rate formulas. There the determining quantity is the product or of 
the cross section with the relative velocity of the reactants. Hence 
the quantity a , which we term the primitive rate coefficient is the 
thing for which we desire an effective mathematical representation. 

In this connection it is noted that Axford, et (1963) suggested 

taking a constant, a choice that can lead to errors somewhat exceeding 
a factor of two when applied over the entire solar wind regime. But 


because the cross section is slowly varying with velocity, roughly like 
Y over the solar wind regime, recent authors (e.g«, Semar, 1968; 

Holzer, 1972; Fahr, 1972) have taken the charge exchange cross section 
to be the constant appropriate to the solar wind velocity at earth. 

The latter choice of itself, implies errors of less than a factor of 
two. But taking constant implies the mathematical problem of 

treating the relative velocity, in the past with some approximation that 
introduces additional error into a . Whereas taking o' constant 
obviates the necessity to prescribe the relative velocity. Using equa- 
tion (2) we investigate analytically the question of making an effective 


formulation for a • 
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First it is noted that the coefficient of the last term in equa- 
tion (2) is small offering the opportunity for effective linearization 
in V . To this end we consider variations in velocity about 

some constant velocity . Then equation (2) may be expressed as 


= 62.3 - 33.5 logj^Q + 4.5 log^^ (1.4b-3) 

- (33.5 - 9.0 logj^Q v^) logj^g ) + 4.5 logJgCj ) . 

, o o 


Taking = 40 corresponding to the velocity of the supersonic solar 
wind at earth and expressing the result in natural logarithms we have 


from equation (3) 

2 

= 20 - 8.25 >fn ^ + 0.85 ^ (1.4b-4) 

e e 

In the subsonic solar wind we base the cross section on = 10 
which is characteristic of both the post shock velocity and thermal 
speed. The choice yields 


4v 2 4vv 

= 33.3 - 10.6 (^) + 0.85 m (— ) 


(1.4b-5) 


Lastly in the interstellar wind regime we base the cross section on the 
characteristic thermal speed for which v^ = 2 , In this case we get 


= 52.6 - 13,4 fjn 






'/XiU 


(1.4b-6) 



•-() n| CO I 
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Equation (4) , (5) , and (6) have the form 


2 V 


a - Oq “ "^1 ^ V + ^2 ^ V 
o o 


(1.4b-7) 


Expanding in Taylor series about gives, xHth Av - v - , 


o = a - 


= o - 


a, . 2 

Av , ^ , V . Av , 1 / 


(1.4b-8) 


* (Oi + 4a«) » 2 

Av . M 2^ ^Av^ , 

O- — + o ( — ) + 

1 V 2 ' v 


Finally multiplying equation (8) by v gives the series for the 
primitive rate coefficient 


Av 


(£>1 + 2 

or = (<j^ - ^ + ''o 2 ‘ 

0 o 


(1.4b-9) 


In the equation the symbol represents the. product v^ cr^ • 

Table 5, Parameters of formulae for primitive rate coefficients 
for charge exchange 


Region 

1. 

V 

o 

2. 

a 

o 


d 

(p “cr. )v 

■ O 1“^ 0 

(p^+4a2) 



supersonic 
solar wind 

40 

20.0 

8.25 

8.04 

4.70 

2.33 

subsonic 
solar wind 

10 

33.3 

10.6 

3.33 

2.27 

0.70 

interstellar 

wind 

2 

52.6 

13.4 

1.05 

0.784 

0.17 

~ velocity in units 

— cross sections in 

— products of cros s 

g 

of 10 cm/s 2 

units of 10 cm 

section with velocity 

have units of 10 

3 -1 
cm s 


V 

o 
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Table 5, is a compendium of parameters of the primitive rate 

coefficient for the three regions of flow. We note first that the. 

terms in the last three columns decrease monotonically to the right. 

Then inspection of equation (9) reveals that in any sub region where 

^ « 1 the first term ^ in the series provides an excellent approxi 
V o 

o 

mat ion to the primitive rate coefficient. But more to the point, in any 
A\f 

region where — « 1 , a good approximation is provided by the linear 

o 

representation involving the first two terms. And in the context of the 
overall development of the subject, taking only the first term could be 
considered adequate. In the succeeding subsection c. we shall make 
use of the linear representation for o? in a linear thoery for the 
supersonic solar wind region. 

The reaction rates per particle for protons and hydrogen atoms are 
given respectively by 

or'^=aN and (1.4b-10) 

c 

Of = Of n (1.4b-ll) 

c 

where for protons N is the number density of target hydrogen atoms 
and for hydrogen atoms n is the number density of target protons. 

The corresponding mean free times against charge exchange for a proton 

■f* 

or a hydrogen atom are given respectively by the reciprocals of (X^ 
and a • Finally the volume rates for charge exchange reactions for 
protons and hydrogen obey the symmetric relations 

n - n=OfNn and (1.4b-l2) 

c c 

ft 

N = aN=anN 
c c 


(1.4b-13) 
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For charge exchange in the interstellar wind region, the flow is 

6 * ■ 

subsonic and the characteristic speed 2 x 10 cm/s is the square root 

4 o 

of two times the mean square velocity for hydrogen at 10 K . From 

-8 

Table 5. the primitive rate coefficient *5^ has the value 1 x 10 
cm'^/s . From equation (10) the reaction rate and mean free time for 
protons against charge exchange in the interstellar wind are then 

1 X and 1 X 10^ s respectively; and for hydrogen atoms, 

^ ,^-10 -1 , - -„10 

1 X 10 s and 1 x 10 s . 

Based on the mean free times we can construct some revealing mean 
free paths for interactions. First for the protons in the ionozed 
interstellar* wind, the characteristic speed for bulk transport between 
charge exchange interactions is the bulk velocity 6 x lO^cm/s , The 
mean free path for bulk transport of protons is given by the pro- 

duct of the bulk velocity with the mean free time. The result is 

+ 14 , 

X - 6 x 10 cm or 40 a, u. 

The cross section for hydrogen atoms to scatter hydrogen atoms is 

2 

given approximately by ^ira^ where a^ is the first Bohr radius 

•*8 “16 2 
,529 X 10~°cm . The cross section is thus roughly 3,5 x 10 cm , 

Taking the reciprocal of the product of the cross section with the 

-3 

hydrogen density 0,1. cm provides the self scattering mean free 

path « 3 X lO^^cm or 2000 a.u. 

Lastly taking the product of the characteristic thermal speed 

6 

1,5 X 10 cm/s with the mean free time for charge exchange gives the 

mean free path X^ for hydrogen to interact with the ionized inter- 

16 

stellar wind as 1,5 X 10 cm' or 1000 a.u. 
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The characteristic dimension of the flowfield for the interstellar 
wind can be taken as the distance from the sun to the heliopause boun- 
dary, roughly 300 a.u. Comparing this value with the mean free paths 
for protons and hydrogen atoms we observe a striking dichotomy. Speci- 
fically, as was stated in section 2, , the mean free path for hydrogen 
to collide with either itself or the protons of the interstellar wind is 
long compared to the dimensions of the flowfield. This result permits 
the use of collisionless kinetic theory to describe the transport of 
interstellar hydrogen into the solar wind. 

On the other hand, the mean free path for charge exchange in the 
bulk transport of protons of the ionized interstellar wind is small 
compared to the characteristic dimension of the flowfield. Consequent- 
ly the protons of the ionized interstellar wind frequently will be ex- 
changed for new ones obtained from the neutral wind. But the ionized 
interstellar wind, flowing as a continuum fluid around the heliopause, 
does not have the same mean velocity as does the neutral hydrogen which 
is not significantly deflected by the flow of the ionized component. 
Hence each charge exchange event in the interstellar wind in the vici- 
nity of the sun results in an exchange of momentum contributing to a net 
impulse on the ionized interstellar wind. Quite evidently the impulse 
is in the direction of the bulk velocity of the neutral hydrogen and 
results in pressure of the ionized interstellar wind on the heliopause. 

The Mach number of the interstellar wind is low and the random 
velocities of the protons exceed their bulk speed. Consequently what 
small percentage of the primary neutral hydrogen does suffer charge 
exchange close to the sun will to a large extent be replaced by sec- 
ondary hydrogen with nearly the same velocity distribution as that of 
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the primary hydrogen. Hence the flux of hydrogen across the helio- 
pause into the subsonic solar wind is hardly altered by the presence 
of the ionized interstellar wind. Thus in evaluating the transport 
of neutral hydrogen into the solar wind (Chapter II.) we assume at 
the heliopause boundary the components of the differential flux of 
hydrogen directed into the solar wind region to be given by the 
Maxwellian distribution with temperature and density characteristic 
of the undistributed interstellar medium. 

For charge exchange in the subsonic solar wind the characteristic 
speed of relative motion is the mean velocity 10 cm/s for the tern- 
perature 10 of the subsonic solar wind. From Table 5. the pri- 

-8 3 “1 

mitive rate coefficient is 3.33 X.IO cm s , Following equations (10) 

and (11) with densities from Table 4., we find the charge exchange rate 

-9-1 8 

and mean free time for protons to be 3.33 X 10 s and 3 X 10 s ; 

-II -1 10 

for hydrogen, 3.33 x 10 s and 3 x 10 s , 

Looking first at the numbers for hydrogen we see that the mean 
free time 3 x lO^^s is the same while the characteristic speed is 
thrice as great in the subsonic solar wind as in the ionized inter- 
stellar wind. Thus the mean free path ( ^ 3000 a.u.) is three times 
long in the subsonic solar wind as in the interstellar medium. 

Hence the presence of the subsonic solar wind also does not greatly 
effect the transport of interstellar neutral hydrogen into the super- 
sonic solar wind, j 

Next considering the mean free time for subsonic solar wind pro- 

^15 

tons, we have that the corresponding mean free path is 3 x 10 cm 

7 

based on the post shock bulk speed of 10 cm/s , Since the mean free 
path ( « 200 a.u.) is comparable with the characteristic dimension 
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of tho subsonic region, as measured by the thickness ( 180 a,u.) for 

example, each solar wind proton has a significant probability of suf- 
fering a charge exchange event in the region. 

For charge exchange in the supersonic solar wind the character- 

istic speed is the bulk velocity 4 x 10 cxn/s \ From Table 5, the 

-8 3 

primitive rate coefficient is 8 X 10 cm /s , From equation (10) 

with Table 4., the charge exchange rate and mean free time for protons 
— q 8 

are 8 ,X 10”^s and 1,2 x 10 s respectively. Similarly for hydrogen 

* _2 -1 

the rate per atom is 4 x 10 (r ) s and the mean free time is 
6 *2 

2 X 10 r s . 

Based on the mean thermal speed 1.5 x 10 cm/s , the mean free 

12 “^2 ^2 

path for hydrogen is 3 X 10 r cm or 0.2 r a.u. Virtually all 
of the hydrogen that passes within a few a.u. of the sun will suffer 
charge exchange while beyond ten a.u. hydrogen is not greatly affected. 

Based on the bulk speed, the mean free path for protons is 5 x lO^^cm 
or 300 a.u. The radius of the supersonic region being half the mean free 
path, a significant fraction of the solar wind protons will undergo charge 
exchange in the supersonic region. 


\ ■ . ' • 

Table 6. Mean free paths for charge exchange and characteristic 

flowfield dimensions by region. 


supersonic 
solar wind x_ 

subsonic 
solar wind 

interstellar 

wind 

proton mean 
free path 

. 300 

200 

40 

hydrogen mean 
free path 

0.2r^ 

3000 

1000 

characteristic 

dimension 

145 

185 

300 

lengths in astronomical 

units 



A convenient comparison of mean free paths with characteristic lengths 
for the relevant regions can be made in Table 6, whose contents summarize 
the previous discussion. We now address the physical consequences of 
charge exchange in the solar wind regions. 

Since (1) the mass of the electron is small, (2) the square root of 

m 

the cross section for the reaction is large 5A) compared to the size 
of the bound atom thereby providing for large separations of the nuclei 
in passage, (3) the electron acts to electrostatically shield the nuclei, 
and (4) the kinetic energy of relative motion is large compared to the 
potential energy of the electron in the neutral atom, the exchange of 
momentum between nuclei in the charge exchange encounter is very slight. 
Thus while exchanging the electron, the nuclei may be taken to change 
neither speed nor direction as a result of the encounter. 

In the act of shocking the solar wind, as vze have seen, exchanges 

momentum flux for pressure, a process which on the particle level serves 

to randomize direction of motion but leaves substantially unchanged the 

characteristic high velocities ( .. lO^cm/s) of the protons. But in 

the same heliocentric reference frame, the characteristic speed of the 

6 

interstellar neutral hydrogen is only 10 cm/s . Hence the result of a 
\ 

charge exchange event in, either the supersonic or subsonic solar wind 
regions is the exchange of a high speed proton for a low speed one and 
the net loss of a lot of particle kinetic energy from the solar wind. 

The loss of energy per event of course resides in the secondary or 
"fast** hydrogen atom. 
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If the event takes place in the subsonic region, the loss in proton 
kinetic energy is primarily reflected in a reductiort of bulk internal 
energy and temperature as the result of the random motion there* The 
pressure being roughly a conserved quantity in the subsonic solar wind, 
the loss of internal energy is balanced there by gains in density and 
magnetic field strength. 

If the event takes place in the supersonic region, the net loss 
in particle momentum is reflected in a corresponding loss in bulk mo- 
raentum flux due to the highly ordered nature of particle motion in the 
hypersonic radial flow. In the light of the shock termination model 
of section 3,, the loss in momentum flux is seen to result in a re- 
duction of the radial distance to the shock as previously estimated 
in the absence of this additional factor in the pressure - momentum 
flux balance. We estimate the effect and others due both to photo- 


ionization and charge exchange In the following. 


c. Linear theory for the supersonic region 

With the omission of magnetic stress -and hydrodynamic pressure 
terms which together only amount to about 3% of the momentum flux of 
the supersonic wind at earth, the momentum equation in conservation 
law form for steady flow with sources (Chapter II) may be written 


div (pv v) - n p 


(1.4C-1) 


where p is the net momentum gain to the solar wind per charge exchange 

• ‘ / 

event. If we neglect the 17. momentum contribution of the primary hydro- 
gen atom and assume spherically symmetric radial flow then from the pre- 
vious discussion p - - m v e and equation (1) may be expressed 
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^2 m V « (X-^c-2) 

r 

Substituting equation (l,4b-12) for the volume rate of charge exchange 

2 

and making the identification p v r - m cp for the mass flux permits 
equation (2) to be cast 


+ . (1.4C-3) 

dr cp dr 


From the results obtained in section 4a. (equation (8) we express the 
proton flux in equation (3) by the approximate linear relation 


9 = <Pe ( + ^ ) = <Pe 

9e 


(1 + 


Of N 
e e 


n 


(1.4c -4) 


Evidently in the bracketed terra of the linear expression we have 

No 

dropped a quantity ^ — . r which is small ( ^ 0.001) compared 

n V e 
_ e e 

to unity. For a , use is made of the linear approximation 


Of = c 7 V + (a - a. ) (v - v ) - o, v + (a - cr. ) v 
oo o 1 o lo o 1 


(1.4C-5) 


where the parameters are to be chosen from Table 5. for the supersonic 
region. Then equation (3) becomes 


dv 

dr 






N + 


!S£1 

n V 
e e 


(1 + 


a N 
ep o 

n V 
e e 


-1 


r) 


V + 




= 0 (1.4c-6) 


Lastly we take N = which is a good approximation over most of the 

supersonic region and treat equation (6) as a first order perturbation 

"ep^o 

problem in the small quantities — r , (a^ - r , etc, 

“e e 
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^ep^o 

Accordingly, we expand the binomial expression (1 + — r) in 

e e 

power series to first order and are left to solve an ‘equation of the 
form 


dv - 2 

^ r) V = 


- a^ V 

3 e 


(1.4c~7) 


a N 

where a = (a - 0 ^ N + . 

e e 


or N 
_ ep o 


- — t — and - cr,N 

n V * 3 1 o 


e e 


The solution of the zero order problem is 


-a-r , 

v-c e i+c^ • 
o 1 


(1.4c-8) 


wnen 


the constant c is replaced by the vaiiable fuucLioii A(i.) and 


the solution substituted in equation (7) there results . 


(A^ - a 2 r A) e"^l^ + (a^^- a^ r) + a^ = O'. 


(1.4C-9) 


Finally, the quantity a^^r being small in equation (9) we expand A 
(and k*) and the exponential in power series and set the resulting, 
coefficients of powers of r to zero. Through first order the result 
is 


A. + a- c + a, V - 0 and 
1 11 3 e . 




(1.4C-10) 

(1.4C-11) 


0 
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We have above only two equations in three unknowns ,'A^ , and 

c . The needed third relation is provided by the initial condition 


V = V, 


which, consistent with the linearization of the problem, we 


apply at r = 0 giving 


A + c- v ^ 0 
ole 


(1.4C-12) 


The solution of equation (10), (11) > and (12) is 


ai e 


(1.4C-13) 


2 

, . ^3 ^2 

and IT - T 

1 


Thence with expression (8) and the subsequent discussion, the solution 
of equation (7) is 


V = 


I a, a^ 3 2 I 

' “1 a‘ 1 ®1 


a, a, 

(1.4C-14) 


to first order. When the binomial products are expanded the relation 
is seen to reduce to the exceedingly simple form 

V = V [1 ~ (a^ + a^) r] 


= Vg [1 - r] . 

e e 


(l,4c-15) 

/ 


/ 

/ 


The number density may now be obtained by dividing the proton flux, 

2 . 

equation (4), by the product v r yielding 
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n - n 


e r 


1 + a^r 


1 - (a^ + a^) r 


“ % r t ^ ^ ^2 *** ^3^ 

e 


r a N 

H n -— [ 1 + (a N + 2 ° ) r ] 

e r o o n V 

e e e 


(1.4c“16) 


to first order. 

Next multiplying the shock relation, equation (1.2-17) , by r 
provides the equivalent expression 


7 + 3 
2(Yfl) 


M (pv 



(1.4C-17) 


In equation (17) the product cpv obtained from expression (4) and 
(15) is 

<pv = 9e''e [1 + (^2 • ^ 

s m V [1 - a N r] :(1.4c-18) 

e o o 

to first order. > - 

With numbers garnered from the previous discussion we have that 

17 1 ^ -17 -1 

a N = 20 X iO cm and — ^ = 7 . 5 x 10 cm 
o o n V 


These parameters are proportional to the charge exchange and photo- 
ionization rates used by other authors. Because the photo Ionization 
rate is only about 3/8 the rate for charge exchange, it has been 
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presumed Fahr (1972), Holzer (1972) Axford (1972) that photoionization 
is "relatively” unimportant but "qualitatively" has much the same in- 
fluence on the system. Equation (15), with the numbers above, shows 
that the presumption is true for the velocity. Equation (16) shows 
that photoionization and charge exchange are of comparable importance 
in increasing the density, a fact that has been more of less under- 
stood. But equations (17) and (18) show that, in a different way, the 
presumption is not valid for the momentum flux and, therefore, the 
location of the shock transition. For in the latter case the first 
order effects of photoionization in decreasing the velocity are exactly 
cancelled by the increase in mass flux due to the added particles. 

Hence in so far as first order effects of any relative magnitude are 
concerned only the charge exchange process plays a role in determining 
the shock location. 

There is addiLioually something to be learueu iicj.e Ii-um reiur 's 
(1972) analysis viewed in the light of the foregoing. In part, the 
point we wish to make involves the desirability of employing to a con- 
sistent level of approximation the lowest order set of moment equations 
necessary to describe the phenomena in question. 

In his analysis of the modified shock location Fahr assumed a 
priori that photoionization had negligible effect on the problem and 
took the mass flux to be conserved. As we have shoTim this assumption 
leads to no first order error in the computation of the shock location. 
Further Fahr assumed, as we have done, that temperature effects were 
negligible. He then described the variation of the velocity of the 
solar wind through an energy equation with a source term. 
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Though it is not specifically described as such, Fahr*s equation 
is not a total energy equation in which thermal energy effects have 
been neglected; and his source term is not simply the net energy lost 
to the flow in the escape of the primary proton. Rather, his equation 
is the so called "mechanical energy equation" of fluid dynamics and his 
source term is twice the net energy lost to the flow. Fahr argues 
for his source term by adding to the net loss of energy to the flow the 
additional "work" done by the remaining fluid in accelerating the secon- 
dary proton. Fahr’s result is correct but because of the complexity of 
the underlying mechanisms his method of getting the needed source term 
requires insight that as argued, at least, is not fully convincing. 

By cohtrast in the momentum equation formulation the source terra . 
is the net momentum lost to the flow in the escape of the primary pro- 
ton, while momentum is conserved in the interaction of the indistin- 
guishable secondary proton and the remainder of the flow. 

The connection between the two formulations which have yielded 
identical results is very simple. In fluid dynamics the mechanical 
energy equation is obtained from the momentum equation by multiplying 
the latter through by the velocity, a process which gives directly and 
unambiguously the source term used by Fahr, 

Since Fahr's source term, the kinetic energy loss is twice the 
total energy loss of the flow, it is of interest to identify where the 
remaining kinetic energy goes. We recall that after the charge jexchange 

event momentum was conserved in the interaction between the secondary 

/ 

proton and the remainding flow. The process terminates when the secon- 
dary proton reaches the bulk flow speed. But' this scenario describes 
the "perfectly inelastic collision" of elementary mechanics. In such 



68 


collisions, we recall, momentum is conserved but kinetic energy is 


not j the process is dissipative, converting the lost ‘kinetic energy 
to heat. 

Thus without understanding the underlying mechanisms we also 
attribute the added loss of kinetic energy in the charge exchange 
interaction to dissipative heating of the flow. Indeed in a refer- 
ence frame fixed in the flow, in which the initial velocity of the 
secondary proton is -v , Fahr*s ’*work" may be seen as done against 
friction in slowing down the secondary proton to rest. 

Returning to the computation of the reduced shock radius, we 
find that combining equations (17) and (18) results in the quadratic 
equation 


2 

r 


+ 2 


Y + 3 

4(y+ 1) 


mcp V a N 
^e e o o 



, _ nKp V 

Y + 3 '6 e 

2(y + 1) P3 


0 


(l.Ac-19) 


In the absence of the second term, the equation gives the classical 

solution r obtained previously, equation (1.2-18). Hence equation 
so 

(19) may be written 


2 

r 


+ 2 r 


2 

so 


a N 
o o 


r 



^ (1.4C-20) 


When r in equation (32) is normalized by r^^ , the solution obtained 
from the quadratic formula is 


a N r 
o o so 


/ a N r \ 
(_oo^) 


I1/2 


+ 1 


r 


2 


2 


J 


(1.4C-21) 



•cr N r 

The quantity ^ ^ -- has the value 0.145 based on numbers devel- 
oped to this point. Thus the radical may be expanded in a power series 
giving 


R 


R 

r 

so 


1 


a N r 
o o so 

2 


1 { ^o^o'^so \ , 

IV 2 / + 


(1.4C-22) 


Thus the inclusion of particle processes in the model serves to reduce 
the estimate of the shock radius by some 14% to 125 a,u. from the 
previous estimate of 145 a.u. 
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d. Lyman-alpha scattering and radiation pressure on hydrogen 

The chromosphere of the sun emits a strong, broad line of hydrogen 

« 

Lyman-cy radiation centered at 1215.67 A in the ultraviolet. Reported 

rocket measurements (Purcell and Tousey, I 960 ; Bruner and Parker, 1969) 

2 

place the integrated flux in the line at about 6 erg/cm sec and the line 


width at roughly 0.8A , 

Bruner and Rense (1969) have calibrated the flux spectral density 
profile over the line. The line profile features a broad solar absorp- ^ 

tion minimum of some 0,3 A width about the line center. Based on twice 
the most probable thermal speed, the effective doppler width of the absorp- 


tion spectrum for the bulk of the hot interstellar hydrogen is 0,1, A; 

falling well within the self absorption notch in the emission spectrum. 

11 -2 -1 ®-l 

There the solar flux spectral density is 4,2 x 10 photons cm s A . 
The rate g at which an atom of hydrogen scatters solar Lyman-o' 


photons is given ( 2 ,g. Barth, 1969) by 

, „ , X TTe^ 

S = V ST ^ 

e 


(1.4d-l) 


where nF, is the flux spectral density per unit wavelength , X and 
X 

V are the wavelength and frequency of the radiation, and f is the quan- 
tum mechanical oscillator strength 0.416. Using the quoted flux density 

-3 -1 

we find g to have the value 2.28 X 10 s at 1 a.u. 

The hydrogen absorbs momentum from an essentially plane wave pro- 
pagating radially from the sun and scatters (emits) symmetrically as a 
dipole. Hence in each scattering event a hydrogen atom suf f ers ^n impulse 

that is directed radially outward and equal to the momentum of .the photon 

/ 

hv/c , Thus the net force of radiation pressure in the Lyroan-of line is the 
product of the scattering rate by the impulse per event 
hv 

•"a = T 8 


(1.4d-2) 
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With the value found for g , the force of radiation pressure at 

1 a*u. is estimated to be 1,24 x 10 dyne. Solar gravity contributes 
an opposing radial force 0.99 X lO"^^ dyne at 1 a.u. The 25% difference 
between these two numbers falls within both experimental error in a single 
flux measurement (Bruner and Parker, 1969) and measured flux variation 
over the period of a solar rotation (Meir, 1969). The extent to which 
solar Lyman-(y emission varies over a solar cycle is not known but there 
are indications it is within a factor of two and at least comparable with 
the variation over a rotation (see Banks and Kockarts, 1973 for related 
discussion). 

Since interstellar hydrogen moves with a speed only on the order of 

2 a.u, per year with respect to the sun, all solar flux variations in- 
cluding those over a solar cycle will be largely averaged out in the 
transport of hydrogen into the inner solar system. Thus in the absence 
of any more definitive information and for the compelling mathematical 
advantages the assumption affords, we assume the force of radiation pres- 
sure exactly balances solar gravitation in the transport of interstellar 

atomic hydrogen. In this connection we note that the Lyman-o' flux varies 
-2 

as r when the sun is assumed to be a point source scattering in an 
isotropic medium. Hence if we were to assume the forces of radiation 
pressure and gravity balance anywhere, we just as well assume they bal- 
ance everywhere and the hydrogen particle trajectories become straight 
lines into the solar system. 

The atomic absorption line width in scattering is given in terms of 
frequency interval by Stone (1963). Converted to wavelength interval 


by the relation = — Av , the absorption line width is 

c 


. 8n /ne^ \ _ . _ 

e 


-n. 




(1.4d-3) 
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Hence the atomic absorption line width is very much narrower than the 
0.1 A doppler broadened absorption. band generated over the distribution 
of the scattering atoms. 

The effective solar flux for scattering within an atomic absorption 
line is 


And the definition of the atomic cross section a for scattering is 

s 


8 = ^s 


(1.4d-5) 


With the relation X = c/v , equations (1), (3), (4), and (5) can be 
combined to give 


3 \‘ . 

<j = -5 — f . 
S OTT 


(1.4d-6) 


-12 2 

The cross section has the value 7.35 X 10 cm . 

For purposes of estimating the mean free path for scattering Lyman-of 

C 

photons, we assume the hydrogen to be uniformly distributed over the 0.1 A 
doppler broadened absorption band. Then the density of hydrogen 

atoms occupying any absorption line in the band is roughly 




a 


N 


dop 


(1.4d-7) 




Assuming the hydrogen to have the undisturbed interstellar number density 
-0.1 cm the mean free path is found to be of the order 10 a.u. From 
the way in which the doppler bandwidth was defined, the above number is 
appropriate to photons scattered from the bulk of hydrogen under the peak 
of the distribution function. For hydrogen in the velocity tail of the 
distribution, the number density is considerably less and the mean free 
path is commensurate ly greater. 

The volume rate of scattering is the product of the rate g per 
atom and the hydrogen number density N • Assuming the single scattering 
model, we get the contribution to the observed intensity from sources in 
dr at r to be 

dl = ^ dr . (1.4d-9) 


The total intensity is then given by the integral of equation (9) taken 
over the line of sight. Neglecting the effect of aspherically symmetric 
multiple scattering, we assume the solar flux and, hence, g to decrease 

as r ^ and 
1 = 

Equation (10) finds use in Chapter III where we calculate the intensity 
of Lyman-c^ backscatter from the distribution of interplanetary hydrogen 
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that we find. We now turn our attention to completing the mathematical 
model including the problem of the transport of hot interstellar neutral 
hydrogen into the solar wind. 
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Chapter II, A Mathematical Model 

In chapter I, we have seen that the mean free paths for protons to 
interact with hydrogen in either flow regime are short compared to charac- 
teristic flowfield dimensions. This fact and other plasma properties 
involving short characteristic lengths lead us ultimately to treat the 


ionized fraction of the gas through a continuum model. 

On the other hand, except very near the sun, the mean free paths 
for hydrogen to interact with either itself or the ionized fraction in 
either flow regime are very large compared to characteristic flowfield 
dimensions. This fact leads us to treat the transport of hydrogen 
through a single particle description. 

But the state of the two gases are interdependent through the photo- 
ionization and charge exchange processes. The global formulation including 
the interaction between the two gases is best understood if we begin the 
discussion where the two descriptions malce contact, at their fouuuciLioat» 
in the kinetic theory. 


1, Distribution Functions, Boltzmann’s Equations, and the Plasma 
Transport Equations. 

^ The state of the partially ionized gas in either flow regime, solar 
wind or interstellar wind, is specified by the set of distribution func- 
tions f^ , f^ , f^ for the three species comprising the gas— electrons, 
ions (protons), and neutral hydrogen atoms. The distribution functions 
which are the locally averaged particle density functions in the six 
dimensional phase space of position and velocity are solutions of 
Boltzmann's equation 



d 

axp 




C + s . 


( 2 . 1 - 1 ) 
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We emphasize for the sake of the subsequent discussion that for each 
species of particle there is a separate equation of the type described 
by equation (1) . 

In equation (1) the summation convention is invoked for repeated 

subscript indices. The quantities a are the cartesian components of 

P J 

particle acceleration due to macroscopic fields «• gravitation, radiation, 
electric, and magnetic. The symbols ^ and ^ represent functions de- 
scribing discontinuous additions and subtractions of particles from unit 
volume of phase space. 

The symbol C stands for collisions in which particles of the given 
species interact with particles of either the same or different species, 
all particles being changed in velocity but unchanged in kind by the 
interaction. In fact only Coulomb collisions within and between the popu- 
lation of electrons and protons are important in the problem, and for hydro« 
gen Cji effectively zero. 

The symbol S represents particle processes in which changes of 
kind do occur. Here the processes are photoionization and charge exchange. 
The interactions between hydrogen and the ionized species are described 
by the functions These are intimately mathematically 

reflated and provide the basis for the source terms of the continuum 
equations for the plasma. 

The macroscopic quantities number density n , bulk velocity V , 
and total specific kinetic energy e are given by the first three velocity 
moments of the distribution function 

= j f d\ 


n 


( 2 . 1 - 2 ) 
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V = - r V f d\ = <v> 

« n J ^ 


(2.1-3) 


, V . V ^ , 

= - J £ d V = i < y . V > 

n Z Z ^ ^ 


(2.1-4) 


When the distributed velocity v is ^^ritten in terras of the bulk 
velocity by 


V = V + v’ 


(2.1-5) 


then < V* > ”0 and 


e = ^ + *1 < V* . v*> . 


( 2 . 1 - 6 ) 


Equation (6) shows the decomposition of the total specific kinetic energy 


12 

into the sum of the hulk specific kinetic energy V and the Gpcciti^ 


internal energy, the latter to which we give the standard symbol £ . 

The moment equations for a species result from taking successive 
velocity moments of the associated Boltzmann equation. The continuum 
transport equations in conservation law form for the macroscopic quanti- 
ties n , V , e derive from, the first three moment equations which are 


i-'k = jcd\ + jsd^, 


(2.1-7) 


(n < v >) + 1:; (n < v^v >) ~ n < a > 

Qt ^ a ' dx ^ ^ a ' a 

P 


3 3 

+ -rvCdv+rvSdv 

J or ^ a 


( 2 . 1 - 8 ) 
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+~ (5<VV Vq>) = 

at 4 O' a ' dXg '^2 Of or 8 ^ 

" < Ve ^ + X I Va c 

+ r ^ V V s d^v (2.1-9) 

J 2 O' Of 

In equation (8) for the mathematical convenience we have represented 
the vector equation by the set of scalar cartesian component equations 
for or = 1 , 2 , 3 . 

By virtue of the conservation of particles in collisions, by de- 
finition, the first integral on the right hand side of equation (7) 
is identically zero. Then from equations (7) and (9) multiplied by 
the particle mass m and by using equations (3) , (4) , (5), and (6) 
we get the continuum transport equations for a plasma species. 



where the mass density p = mn . 

In accord with the standard definitions from kinetic theory, in equations 

(10) , (11), and (12) the pressure p , deviatoric stress tensor components 

CT , and conduction heat flux vector components q are represented by 
pa P 



(2.1-13) 
(2.1-14) 
(2.1-15) 

With equation (13) and the definition of the specific internal energy, 
equation (6), we get the state equation 

e = 2 E . (2.1-16) 

2 p 

Lastly a temperature can be defined for the species by 

T = I H e . (2.1-17) 



When n , V , and T are sufficiently slowly varying in space 
and time, the distribution function will approximate a Maxwellian as 
quaranteed by the Boltzmann II- Theorem for the relaxation process. 

Under such conditions the distribution function may be represented by 
a first order perturbation expansion the first, zeroth order, term of 
which is the Maxwellian distribution. By the method of Chapman and 
Enskog the Boltzmann equation can be solved approximately for the first 
order perburbation which is linear in the space derivatives of n , V , 
and T . The approach permits the integration of relations (14) and (15) 
to give approximate expressions for the components of the deviatoric 
shear stress tensor and conduction heat flux vector that are linear re- 
spectively in the derivatives of the velocity and the temperature. 
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In the absence o£ a magnetic field the components of shear stress 
and heat flux take the form 




= d ( 


axc^ axp 


2 , 

3 


dv 


^6 \ 
3x / 


and 




(2.1-18) 


(2.1-19) 




in which the single (scalar) coefficients y, and k are termed the 
viscosity and thermal conductivity, respectively. Using the method 
suggested above. Chapman (1954) has obtained coefficients which may be 
expressed 
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( 2 . 1 - 20 ) 


(z. l-^i; 


( 2 . 1 - 22 ) 


(2.1-23) 


where k — [ 2 1/3/ * ^ 

e n ' 

frtf » 1 as* occurs everywhere in our problem. 

2 1/3 
e n 

When the plasma contains a magnetic field such that the gyto radii; 

/ 

/ 

of particle orbits in the field are short compared to the characteristic 
mean free path for Coulomb collisions, a great deal of anisotropy is 
found in the medium. Then the heat flux vector and shear stress tensor 
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are not described by single scalar coefficients but by sets of coeffi- 
cients that can differ, markedly according to direction relative to the 
magnetic field vector. Braginskit (1965) has derived expressions for the 
transport coefficients in the presence of such a magnetic field. 

When written in terms of our definition gf shear stress components 
and the components h^ of the unit vector h in the direction of the 
magnetic field, Braginskii's expressions become 


"oap "lap 


W o = 4(hh„-^«„)hh w 
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(2.1-24) 
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6^ = 6 ^ - h h* 
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is theKronecker delta, and unit antisymetric 

tensor. 

For protons the coefficients are given by 
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= 0.96 n^KT.T. 
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n.KT.T. 
1 it 


XiCxT+ 2.38) 


rt/ « n.KT.T. 

i , >24 2 r> noN ^ t i 

‘‘1 = ("5 ^ 


n.KT.T. 

4 = 2x.(4x^ +2.38) I 

2i 


(2.1-25) 


where t is the mean free time between coulomb collisions for protons 
i 

with protons, Xj, = identically the ratio of gyro radius to 

eB 

mean free path for coulomb collisions and lu = — ^ is the Larmor or 

i 

cyclotron frequency. Then the parameters and are given by 

Xi + 2.33 and 

= 16 xS 16-12 x^ + 2.33 . 

Finally the mean free time is 

3/2 


3V^^1S^(KT.) 

Ti = 4 

8ne 


where, Spitzer (1962), A- 


(2.1-26) 




(^) 

e 


3/2 


is the ratio of the Debye shielding distance to the impact parameter for 
90° coulomb collisions. The Debye shielding distance 


h - 


1/2 

,KT , 

(2.i 


\4TTn^e / 

/ 


/ 

is the characteristic length over which thermal energy can support local 
fluctuations in electrical charge density in the plasma and is used as 
the effective cutoff distance for coulomb collisions. 


83 


For the electron gas the coefficients of viscosity are given by 


, n KT T 

n^KT^T^ . = (8.2 = 8.50) 


« ,n KT T 2 

= (2.05 Xe + 8.50) , p,® = . (4 x^ + 7.91) ~ 


n KT T 
e e e 


e 7 "e'^'^e'^e 

^4 = - <Xe + 7.91) a 


(2.1-28) 


where very similarly to the case for protons 
+ 13.8 + 11.6 and 

= 18 xf + 55.2 x^ + 11.6 . 


eB 

y ~ w T , u) “ m^c , 
A-e e e e e 


The electron collision time is 


T = 

e 


3j^ (KT ) 

e e 

STTe^n J(/iA 

Q 


3/2 


(2.1-29) 


Similarly, for the heat flux conducted by the electron gas, 
Braginskii’s expressions may be represented 


e eT , ev . 
q = q + q with 


(2.1-30) 


= - xf, h . 7 + X® h X (h X 7 T ) - x * h x 7 T 

q®^ = {P|| h . (v® - v^) - 3 h x[h x (v^-v^,) ] 


G i 

+ 3a^X(v “V)3 n^KT 


e e 
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In the above expressions, the coefficients of thermal conductivity are 




<= (4.66 Xe+ 11-9) 

e e 


(2.1-31 


« n K T T 

e 5 2 « e e e 

- X (o’ X + 21.7) —7^ 

A ■^e 2 A i 


_ m 

e c 


and the coefficients for the additional transport due to collisions 
between electrons and ions are 


P||=0.7 , = (5.1 xf + 2-68) 

(2.1-32) 

3a= Xe(| Xe + 3-05)#^ • 

A'e = x^ + 14.8 Xe + 3.77 . 

\ For the proton gas Bragtnskii gives 

q^. = - H I, h . V T . + h X (h X V T . ) 

- h X 7 T, (2,1-33) 

A ^ .w 1 


) 
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where 


i , 

K,i =3.9 — 

II “i 


n.K T.T. 
1 1 r 


\ + 2.65) 


i /5 2 , , _ 

’‘a = Xi<2 Xi + 2^.65) 


n.K T.T. 
i 11 


..m 

1 i 


and -L' = + 2.7 x^ 0*68 

1 ^1 '^i 


(2.1-34) 


The exposition of terms in the plasma transport equations is 
continued now with considerations of the body force F = m < a > which 
appears in equations (11) and (12) i For the protons the force is given 
by 4 


F = e E + e 


— X B - GM 
c ^ o i_3 


(2.1-35) 


where E and B are the macroscopic electric field and magnetic in- 
duction respectively and e is the charge of the electron and c the 
speed of light, Gaussian units being employed. Similarly, for the 
electron gas the force is 


F ==-eE-e- xB-GMm 
•V ^ . o A. o e ^3 


(2.1-36) 


From equations (7) , (8) and (9) the source terras of the plasma 
transport equations are defined by 


/ 


m 


= Jra s d V 


(2.1-37) 
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P “ r ra V C d^v 


(2.1-38) 


P 

«s 




tn V 


3 

S d V 


(2.1-39) 


®c " J" I v^c d\ (2.1-40) 

®s ° J'l''^Sd\ . (2.1-41) 

o 

The source term m represents the time rate of adding new mass 

bearing particles to unit volume of the plasma. In charge exchange 

the number of plasma particles is conserved and only photoionization 

contributes to the source term. Since there is one proton and one 

electron contributed per event the source term is given by the product 
0 • 

m n where n is the photoionizations rate for. hydrogen (section 1.4a) 



(2.1-42) 


The cource term P^ is the time rate of change of momentum per 
unit volume due to coulomb collisions. Since total momentum is con- 
served in such collisions, there is no contribution to the source term 
from collisions between members of the given species. For the same 
reason, the sum of the source terms for ions and electrons is identi- 
cally zero or 


= - P® (2.1-43) 

Representing the transfer of momentum between the electron and ion 
gases, the source terra gives rise to the electrical resistivity of the 
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plasma in Che one-flu.id description. Where common parameters are the 
same as in the similar representation for the heat flux, Braginski finds 


. . -f 

P® = - b|, h . (V - V*-) - p|, n K (h . OT ) (2.1-44) 

...C II ^ ^ H e ^ e 

+ b h X [h X (V® - V^)] + p n K [h X (h X V T )] 

Xn. A. A mm 

- b^ h X (V* - y^) - n^K(h X 7 Tg) 


In the relations the coefficients of dynamical friction are 


bn = 0.51 


m n 
e e 


(2.1-45) 


(a; - 6.4 + 1.84) ^ 

e e 


' = X.(l-7 xf + 0-78) 


e e 


The source term P , equation (39) , is the time rate of change 
of momentum per unit volume due to photoionization and, in the case of 
protons only, charge exchange. The source term for electrons, due to 
their small mass relative to the protons, does not contribute signifi- 
cantly to the formulation and is omitted here. The source term for • 
protons is given by . 
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== r m.v(o: + O' ) d^v - f ra.vc^^ f.d^v . ' (2,1-46) 

s J p c' N J c X 

In the equation the first integral evidently represents the addition 

of momenta of the new protons arising from photoionization and charge 

exchange with hydrogen. The second integral represents the loss of 

i 

momenta of solar x^ind protons which undergo charge exchange with hydro- 
gen, Retaining only the first term in the expansion for a , we get 
the useful approximate relation (see equations (l,4b-9., 10, and 11)) • 


i rt — 3 “ i 

' “ m V (o' + Of n,)f ^d v - o Np^v 

s p ox^N o 


(2.1-47) 


The source term defined by equation (37) is the time rate of 

change of total kinetic energy per unit volume due to coulomb collisions. 
Since total kinetic energy is conserved among the particles in such colli- 
sions, there is no ccntributior. to the source term from collisiwus belweeu 
members of the given species. For the same reason, the sum of the source 
terms for ions and electrons is identically zero or 



(2.1-48) 


Representing the exchanges of energy between the electron and proton 
gases, the source terms play no explicit role in the one-fluid formulation 
and will not be given here. However the processes which the terms repre- 
sent do play an important implicit role in keeping the temperature of the 
two gases comparable in circumstances where the one -fluid formulation is 


most effective. 



Lastly the source term e is the time rate of change of total 

s 

kinetic energy per unit volume due to photoionization and charge exchange ♦ 
For the electron gas photoionization contributes the excess energy of 
the ionizing photons. As we found in section 1.4a, the average excess 
energy per event amounts to roughly "3 ^ 4.5 e, v, and the source 

term is given by 

- 4 hv n. = i hv cr N . (2,1-49) 

s 3 o 1 3 o p 


Similarly to the case for the corresponding momentum source terra, the 
energy source terra for the proton gas is evidently 


'i r i 2, . P r i 2 +. .3 


in. 


^ r 12, “ ^ t3 — ,, > 1 ..2 • 3 * I 

— I -^v (of + a n.) V - a wn v -r - 5 ^ K T ) (z.j.-oO) 

J2 p oi N ol 2i I i 


Having laid much of the groundwork for the global formulation, we 
now direct our attention in the next section to the one-fluid or mag- 
net ohydrodynaraic formulation of the plasma flow problem. Following that, 
a formulation of the hydrogen transport problem as an approximate solu- 
tion of the Boltzmann equation concludes the chapter. 
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2. The Magnet ohydrodynamic Description 

The magnet ohydrodynamic description of the plasma flow problem is 
a synthesis of Maxwell's equations, the plasma transport equations, 
and equations of state into a reduced set of equations for the state of 
a single ionized fluid and the embedded magnetic field, 
a. Max^> 7 ell's equations and the MHD approximation 
The presence of the rarefied hydrogen gas contributing no signi- 
ficant effect, the electric and magnetic fields which appear in the 
body force terms, equations (2.1-35 and 36), of the plasma transport 
equations are governed by Maxwell's equations for charges in a vacuum 

7 X B 


1 

c Bt 


, 4tt . 

+ — * j 

c t 


(2.2a-l) 



(2.2a-2) 


V X E 


i — 

c at 


(2,2a-3) 


V . E = 4rr(n^ - n^) 


(2.2a-4) 


In the equations, written for Gaussian units, j is the current density 
given by 


3 = e(n^ . 


(2.2a-5) 


The principal assumption of magnetohydrodynamics is that the plasma 
is electrically quasi neutral over the volume element of mathematical 
interest, i.e. one whose dimension is large compared to the scale of 
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charge density fluctuations as measured by the Debye shielding distance 
but small compared to characteristic flowfield dimensions. Debye lengths 
being short, on the order of kilometers, the assumption is applicable in 
our problem and we take the state variable for particle density 

J 



(2.2a-6) 


Then from equation (5) we understand the density of electrical current 
flowing in the plasma to be given approximately by 

J = e n(V^ - V®) (2,2a-7) 

We now and hencforth explicitely assume a steady flow problem in 
the heliocentric reference frame and ignore the unsteady terms in 
Maxwell’s equations and the plasma transport equations. Accordingly, 
the latter combined with the state equation we now write in vector 
form explicitely for each species 

(2.2a-8) 
(2.-2a-9) 


(2.2a-10) 


div (p^y ) - 


„ n.e . 

*dtv = n^eE + -^(V X B) 


—V r + 

_3 ^ ^s 


.1^1 . 5 Pi, - ..i 1- - -1 


div [p.V (-2+ 2 


) - cr_,. V + q ] = n eE . V 


i ^ 


GM p , . • .f • 4 
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div (p V ) = m n 
* e e 


(2.2a-ll) 


n e 

div (p V®V^ + pI-a) = -neE- -^(V® x B) 


(2.2a-12) 


GM p • • 

— ^ r + P + 


div [p^f ( ^ + I ^ . v" + q"] = - n^eE . f 


{2.2a-13) 


GM p 
o*^e 


V^ + 


+ e 


e 

8 


The one -fluid hydrodynamic equations are obtained by adding the 
corresponding plasma transport equations for electrons and protons. 
Thus the continuity equation is 


div (pV) = mn 


(2. 2a- 14) 


where p = "** Pe » ® ®e ’ bulk velocity of the plasma 

is defined by v * i (p.V^ + p V^) , Since m « m. and n — n. « 

p « p and for all intents and purposes p - mn = p and V = , 

e 1 _ !• M M 

With the above definitions, approximations, and assumptions, the 
one-fluid momentum equation is effectively 

o 

M . GM p • 

div (pW +pI-CT)=-jXB ^ r + P (2.2a-15) 

r 

where p=p +p^ ,a=a +cr. » and P-P + p +p +p*‘ . 

e *^i e i c ^,8 c 
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In equation (15) Che omission of a body force term involving the electric 

2 

field introduces an error only of order (V/c) • Again since 

and the coefficients of viscosity (section 2,1) are proportional to the 

masses of the plasma species, it is evident that only the protons con- 

' i • * A# ^ 

tribute significantly to the shear stress and a = a . • From equations 

J " * i 

(2.1 - 43 and 46) and the accompanying discussion we have also £ “ 

Similarly the one-fluid energy equation is obtained as 

div [pV (-pr + “T ^) - a . V + q) == E . j - — ^ r + e (2.2a-16) 

G i * ^ c * i 

where q = 4* q and e = e + e (see equation 2.1-48.). 

^ ^ s s 

When it is assumed finally that the plasma is in local thermodynamic 
equilibrium at temperature T = , there results the state equation 

p = 2 uKT ■= 2pftl (2,2a-'17) 

where ft is the universal gas constant. 

The set of six equations (1) , (3), (14), (15), (16), and (17) are 
now found to involve the seven unknowns E,B,j,p, V,p> and T . 
The equation needed to form a closed set is obtained by combining equations 
(6) , (7) , and (17) with the equations of motion for the electrons equa- 
tion (12), As the result of the small mass of the electron, the viscous 
and inertia terms, except the source term involving dynamic friction 
from Coulomb collisions, may be dropped in the latter equation. Since 

in our problem v = u) T is alv/ays several orders of magnitude greater 
'^e e e 

than unity, the source term for dynamic friction, equation (2.1-44), 
greatly simplifies to 
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m -n . . 

« (v^ - V^) (2.2a-18) 

m n . 

+ 0.49 h . (V® - 

T ~ -V 

6 

- 0,7 n K h . 7 T 

e e _ # 

Then the final result we seek may be expressed as a relation for the 
electric field 

E = - ivxB + — j X B - ~ V p + T] j (2,2a-19) 

^ c ^ nec i ^ 2ne ^ t 

- 0.49 ^ j , BB - 0.7 (Vr) . BB . 

B " eB ~ 

Equation (19) is regarded as the generalized Ohm's law for the plasma; 
the electrical resistivity T[ appearing therein is given by 

m 

11= -4— . (2.2a-20) 

ne Tg 

The terms on the right hand side of equation (19) are ordered from 

the left roughly according to importance, as determined from dimensional 

♦ 

analysis based on characteristic parameters of the flowfield from Chapter I. 
Except in the immediate vicinity of the stagnation point where the velocity 
is very small, the first term is orders of magnitude larger than the others. 
The second and third terms are of comparable importance followed with lesser 
importance by the remainder. Then as a matter of mathematical convenience 
that introduces no inordinate error we drop the last two small terms that 



involve anisotropy. 



95 


For steady flow, the first of Maxwell's equations, equation (1), 
implies 


div j = 0 . 


(2.2a-21) 


Then substituting equation (19) in (3) and using equations (1) , (2) 
and (21) we get 


t(B .V)v X B - (7 X B) . VB 


1 1 
+ - 7 p X 7 -T - - V p X(B. 7B) ] 
P /V M ^ P A/ M n/M 


(2.2a-22) 


- (V.y)B -(V.V)B +(B.V) V = 0 

The last three teruis iu equatloti (22) derive from the generally dominant 
first term of equation (19) and taken alone describe the apparent trans- 
port of the magnetic field lines with the fluid in the so called frozen 
field approximation. Where the approximation is valid the equation may 
be written succinctly 

\ •. . • ■ . 

a 

V X (V X B) = 0 

Finally, from Maxwells equations the vector identities for steady flow 


h (!! -T^) 


E . j = - ^ div (E X B) 


and 
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pemiit the momentum and energy equations, equations (15) and (16), 
to be cast 


_2 ^ BB 

div [PW + (p + 1^) I - ~ - oI = 


GM p 
o 


r + p 


(2.2a-24) 


div [pv + I £) + E X B - V . CT + q] 

GM p 
o 

“ ‘ 3 

r 


(2.2a-25) 

V + e 


Equations (1) and (19) can be used to eliminate j and E from 
the momentum and energy equations. Then collecting from section 2.1 
transport and source terms appropriate to strong magnetic fields » 1) , 
we have the hydrodynamic equations. 


div (pV) = mn 


2 

. r 

n = OTpe ^ N 


,2 ^ 


BB 


div [pW + (p + ^) I - or - = - 


4n 


GM p 
o 


3 :: 


r + P 


a = a. . € .e , 
i-J j 


ij ^o oij 


W . . - I b'^(b.b. 

oij 2 ^ r j 


- I 6. .) B B W . i 
3 ±y a ^ ap 


dV bV 2 

~ ^ ^ ** 6 — 

ij dx. 3x 3 ij dx. 


pKTr 

M a 

^o ra 


(2.2a-26) 

<2.2a-26a) 

(2.2a-27) 

(2.2a-27a) 


/ 
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eB 

X - — T 
^ me 




#r = -..v .!■■ 1.I ■ ■ rf~ ' A _ 

Stt^Ap ’ 2Vm/ 


2 - 
r p 


P = r m v(Q'__ — ^ H f-,d V - O' NpV 

^ J pe 2 m ^ N- o ^ 


(2.2a-27b) 


div [pv(^ + I f •'■ li;p> - • <<’ + i 


4tt 


BB) 


(2.2a-28) 


+ «l;;z B X V (p + b£) 1— B X(B.7) B 

8nne . . (4n) » 

,.2 2_ GM p 

• ! to fto)2 ® • ! 5 + 3^ = ■ l-l + ® 


q = - -| BB . vr - 1^2 BB . (V X B) KT 

^ HlitTC <s/^ ^ ^ 


(2.2a-28a) 


./ = B 9 


PKTt 


3/2 1/2 

m m 

e 


e == 


2 - 
r O' 
e . o 


nm 2 *e,^o v£j 3 

1 2 < V 2 + p) V '' 


(2.2a-28b) 


- O' Np i hv a N 

o 2 2m ' 3 o pe ^2 


When the distribution function for hydrogen is regarded as known, 

equations (22) or (23), (26), (27), and (28) with the equation of state, 
equation (17), form a complete set in the unknowns B , p , V , p , and 
T * We refer to the set as the magnetohydrodynamic equations. 


b. Discontinuity relations of the MIIP formulation 
The Maxwell's equations (2) and (3) and the hydrodynamic equations 
(14), (24), and (25) admit respectively the following integral relations 
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across a surface of discontinuity 


'Vi - ° 


(2*2b-l) 


[E X n] = 0 
-'I 


(2.2b-2) 


[pV ] = 0 

1 


2 

[pw^ + (P + ^) n - ^ 


4tt 


cr , n]^ = 0 


(2. 2b -3) 


(2.2b-4) 


I In ® ° (2-2b-5) 


In the relations, n is a unit vector normal to the surface of dis- 


continuity and B = B , n , etc. 

n ^ ~ 


V’hen only the first term is retained in eejuation (19), i.e* 


E = - - V X B 

c ^ * 


(2.2b-6) 


the substitution made in equations (2) and (5), and the viscous terms 

in equations (4) and (5) dropped, there results the set of approximate 

# 

relations 


[B ]" = 0 
n 1 


(2.2b-7) 


[B V - V B] = 0 
‘ n^ .n.,, ^ 


(2/2b-8) 


fpvj^ = 0 


(2.2b-9) 
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[pV^V + (p + 


8n 


) n - 


B B 2 
.Jlrl 
4n^l “ 


[PV^(^ 


5 p 

+ ^ £ + 
2 p 


|- ) - i B B . Vl^ 
4ttp 4tt I 


== 0 


(2.2b-10) 

(2.2b-ll) 


The normal component of equation (8) vanishes identically leaving for 
the component tangential to the surface 


[B - 
n^t 




(2.2b-12) 


For the normal component of equation (10) we have 



+, p + 



n ,2 

8n a 


0 


(2.2b-13) 


while for the tangential component. 


(pV - X B B if = 0 
^ n t 4 tt n L 1 


(2.2b-l4) 


Two classes of solutions of equations (7), (9), (11), (12), (13) , 
and (14) are of interest. The tangential discontinuity relations and 
the quasi-normal shock relations. The former class of solutions is able 
to describe conditions across the surface of contact between flows of in- 
dependent character and origin. Evidently in our problem such relations 
should be applied at the boundary between the solar and interstellar winds. 
The latter class of solutions apply tooonditions across the shock in the 
supersonic solar wind. 

For the tangential discontinuity relations, as a. matter of definition 
V 

n 


= 0 . 


(2.2b-15) 
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Then from equations (7), (12) and (14) we demand 

B ^0 (2.2b-16) 

n 

since the alternative demands (of the "contact" discontinuity relations) 
[Vj.]^ = 0 and [B^l^ = 0 


would be too restrictive to describe the phenomena of interest. Equation 
(11) is then satisfied identically by equations (15) and (16). Lastly, 
equation (13) provides a total pressure balance condition 


®t 2 


(2.2b-17) 


Retrospectively and for consistency w5_t*b the th#> ■'Vis- 


cous terras omitted from the analysis should be small and therefore 
approximately satisfy the compatibility conditions 


^nn 1 ^nn2 


°tn + “In’? = ° 

or they should be retained with the other terms. In the light of equa- 
tion (16) and Braginskii^fe finding (section 2,2a) that the viscous shear 
stress is greatly reduced in directions along and across strong magnetic 


(2.2b-18) 

(2.2b-19) 

(2.2b-20) 
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field lines and similarly for the conduction heat flux n-orraal to field 
lines, the terms in equations (19) and (20) are indeed small and their 
omission from the ensuing development is warranted. But if viscous 
terms are to be retained in the overall formulation then equation (18) 
•is not justified and instead equation (17) is replaced by 


[P + 



2 

nn 1 


0 . 


(2.2b-21) 


The relations that we shall now derive and term the quasi-normal 
shock relations hold when the momentum flux greatly exceeds the stresses 

2 j 1- 

of the magnetic field (pV » -^) in the supersonic region and the 
flow approaches the surface of discontinuity at nearly normal incidence. 
Such is believed to be the case in the solar wind. We begin the develop 
ment with relations (7), (9), (12), (14), (13), and (11) which can be 
■Written respectively 


B = const 
n 


pv 


n 


const 


pV [- = B [V 1? 

n p X X n 


B 


■ ,2 




? “t 2 




(2.2b-22) 


(2.2b-23) 


(2.2b-24) 


(2'2b-25) 


/' 


(2.2b-26) 


(2.2b-27) 
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By hypothesis the right Iiand sides of equations (24), (25), and 
(27) are second or higher order small quantities. Accordingly, we seek 
as an approximate solution a solution to the set 


2 2 

[V = pV [i = 0 

n^t^l ^ n p ^t 1 

2 2 

[pV + P + ■?“ ], = 0 
n 8n I 






0 


J 


(2.2b-28) 


(2.2b-29) 


(2.2b-30) 


along with relations (22) and (23), Within the scheme equation (25) 
is now regarded as providing the jump in V after B is found from 
equations (28), (29), and (30) * 

When - Vjj * J - pn , and = uB^, then equations (28), (29) , 
and (?0) may be expressed explicitsly in the form 


= const 


2 2 

P2 = Pi + (1 - » J H- 2 (1^) 

1 P 


^ 

■2 2 J 4njpu 


0 


2 

^*1 5 

— J: + 
2^2 


Pi“i 


_0 

4nJu, 


where the parameter p is defined 




(2.2b-31) 


(2.2b-32) 


<2.2b-33) 


(2. 2b -34) 


Substituting equation (32) in (33) provides the equation 
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+ I S ( 1 - g) +1 o 


1 ) 7 ^ + 

Ju^ 


(2.2b-35) 


7 [|(0^ - 1) + a - g)i 

P 




4tiJ u' 


Equation (35) has two roots 3 1 


3 - g [1+5 )] 


1 8nJu. 


8nJu, 


1/2 


(2.2b-36) 


The former solution represents continuous flow and is of no interest 
here. For the latter solution in the supersonic solar wind region both 

2 

« 1 . 


b; 


'’l - ’’l _ V. 

— = — 2 « 1 and —3 =. ^ 

1 p^Uj^ SttJUj^ 


Then expanding equation (36) in the appropriate convergent power series 
and retaining only terms to first order yields the strong shock solution 



B 


[1 + 5 


+ 9 


tl 


P«i 




] = ~ 
2^ 4 


[1 + ^ + 

M, 


2M 


Atl 


(2,2b-37) 


In the latter representation of the solution and sre respec- 

tively the ordinary preshock Mach number and the preshock Alfv^n Mach 
number based on the transverse component of the magnetic field. 

Thus from equations (34), (22), (23), (28), (25), and (29) with 
3 given by equation (37) we have the complete approximate solution 


n2 


= 3v 


nl 


(2.2b-38) 
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B « = B , 
n2 n L 


(2.2b-39) 


*^2 “ 3 


(2. 2b -40) 


^2= 3"tl 


(2. 2b -41) 


B 


^t2 " ^tl ^ p . ^tl 


(2. 2b -42) 


P2= p^+ (1 - 

P 


(2. 2b -43) 


As was the case for the tangential discontinuity relations, if 
viscous terms are to be retained in the formulation then it should be 
sufficient to add the normal stresses to either side of equation (43). 

It is expected that the omission of shearing stresses is appropriate 
in the region of nearly normal flow where the above relations apply, 
c. Equations of quasi -one -dimensional compressibl e flow 
With certain assumptions it is possible to construct an essentially 
one dimensional flow problem that embodies major features of the termi- 
nation of the solar wind. With certain approximations the resulting 
mathematical formulation involves an equation which in the absence of 
magnetic fields reduces to the one -dimensional compressible flow equation 
of gas -dynamics. 

In the interstellar wind region, the assumptions are that ±ji the 
uniform flow far upstream of the sun the velocity vector is parallel to 
the solar equatorial plane and the magnetic field is orthogonal to the 
velocity. In the solar wind region, symmetry is assumed with respect to 
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the solar equatorial plane. Then the regimes of quasi-one -dimensional 
flow are about the upstream and doimstream axes through the sun. 

The approximations are that E = - — V X B and the viscous terras 
in the hydrodynamic equations can either be ignored or treated as pre- 
scribed functions in the solution of the equations by iterative methods. 
Then the mass equation (2,2a-26) 

div (pV) = mn (2.2c-l) 


with the familiar vector identities permits the momentum and energy 
equations (2.2a-27 and 28) to be cast respectively 


2 2 
p7 I pV X (7 X V) + V (p + Ij-) -- 5 ^ div BB 


« pg + div a + P - mnV 


( 2 , 2 c- 2 ) 




2 

- ^ + Ir.) V . Vp + I pV , V (|) - i, BB : W 


V 4np' : 


4n 


= p g. . V + (div cr) . V + a : W 


♦ # 


+ € + Q " inn (• 


2 2 

" + 2 p + ^p> 


(2.2C-3) 


GM 

where g = - — r and Q = - div q . When equation (2) is dotted 
with the velocity and the result subtracted from equation (3) the 
difference is the generalized heat equation 
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• 2 2 

•I pv . V (-) + V . V. I - ( £ + y „) V . VP - i BB : TV 

2 ^ P - ~ 8n p 4rrp^ ^ 4n 




• ♦ • 


= e 


2 2 

+ Q - P . V + m. (|-|£-|„p) + a : W 


(2.2C-4) 


Finally with equation (2.2a-2), equation (2,2a-23) may be written 


B.W-V.VB -(V.V)B=0 


(2.2C-5) 


For one dimensional flow equations (1), (2), (4), and (5) find 
expression 


lim 1 d ■ . . 

A ^ T •j (puAs) = mn 
As dx ^ 


(2.2C-6) 


du , dp 1 dB^ , 1 B^ 

pur: H-r^+ TT.-T— + T* — 


d(j. 


dx dx 8tt dx 4n r 


B 


gp + P - mnu + 
X 


XX 


dx 


(2.2C-7) 


2 2 2 

3 d -p. ^ u dB /P . B B u 

2 dx ^ 0 ^ 8tt dx ^ 4 ttp^ dx 4 tt r, 


=** € + Q - P u + mn(^ “ 4 ? - (2.2C-8) 


2 2 p 4np^^ ^xxdx 


« u , dB lim Id / A \ B n 
3_ -u^~ (uAs) B = 0 , 

r^ dx As-O As dx 


(2.2c-9) 


In the equation As is a differential element of area normal to the 

axis of flow and r and r are the local radii of curvature of the 

B V 

magnetic field lines and the normal surfaces to the velocity field 
respectively. 
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We approximate both and by ^ . Then 

ceasing the repetitious writing of the limit tiS-O , *we can rewrite 
equations (6) -(9) in the forms 


i l£ = 

p dx 


i ^ 1 d As ^ ^ 

u dx dx pu 


(2.2C-10) 


du 

U — “ 

dx 


1 d p ^ 1 dB"^ . B 

p d^ ^ -^ + 


^ dAs 

8np dx 8rrp As dx 


P « ..da 
X U , 1 XX 

- - g H mn — + — 

P P P 


dx 


(2.2C-11) 


2 d p 1 ^ dAs 

2 dx^p^ 8np dx '‘p 4np p dx " 8np As dx 


(2.2C-12) 




1_ 1 ...— 1/2^ fi 

-T — 1 /9 dx ' ' . 


A51/2 


^ 1 ON 


Equation (13) provides an immediate integral which squared and then 
differentiated yields 


\ l_ 

8tt dx 


B^ 1 du .1 dAs 


4tt u dx 8 tt As dx 


(2.2C-14) 


Differentiating p = (— ) p we have 

‘ P 


1 ie = 1_(P\ + /£■) i i£ 

p dx dx^p' ^P"^ p dx 


(2{2c-15) 


In equations (11) and (12) we substitute for the quantities on the 
left hand sides of equations (10), (14), and (15), In the resulting 
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equations we then replace — by ~ a and 

P ^ 

2 

A are the squares of the speeds of sound 
The results of doing the above are 


B 2 2 

■T — by A where a and 
4t7P 

and Alfveh waves respectively. 


du 3 da^ u ^ ^ d &s 3^ a^ 

dx 5 dx 5 ,2 dx 5 As dx 5 Pu 
M 


u du 

- -2 dJ = - 8 + 


M. 


' . da 

X u ,1 XX 

— - — mn + — 

) p p dx 


9 da ^ j. ^ ^ ^ d As 3 ^ 

10 dx 5 ^,2 dx 5 As dx 5 pu 

M 


pu 


{ e + Q 


• u" 3 2^ ^ 

P^u + mn(-^ - ~ a ) + a 


du>. 

XX dx-^ 


(2.2C-16) 


(2.2C-17) 


Where m = is the Mach number and M. = t is the Alfven Mach number, 
a A A 

2 

Finally when equations (17) is multiplied by -j , subtracted from equa- 

2 

*'■* 0 *^ the resulting difference divided by u , we arrive at the 

equation for quasi -one -dimensional compressible flow 


1 

1 ^ 

1 du _ 

1 

-L 

d As 

g 

2 


^ u dx 


As 

dx 

2 

M 


M 


u 

4 

« 

mn . 

5^x 

2 

• • 

e-K) 

. 1 

da 

. XX 

3 

— + 
pu 

3 2 
pu 

3 

3 

pu 

**■ “ 
pu 

^ dx 


(2.2C-18) 


^ \ 
3 u dx ^ 
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3. Transport of Hydrogen of Interstellar Origin and Moments of the 
Distribution Function 

The transport of interstellar atomic hydrogen into the solar wind is 
governed by Boltzmann’s equation (2.1-1). In the equation the self colli- 
sion term can be neglected as the result of the associated large mean 

J 

free path compared to the characteristic dimension of the solar wind cavity 
(section 1.4b).. In section 1.4d we have provided a rationale for omitting 
from the equation the forcing terra involving macroscopic fields. Finally 
in the heliocentric reference frame in which the problem is regarded as 
steady, the Boltzmann equation for hydrogen is reduced to 


div 





(2.3-1) 


Since the fast secondary hydrogen that is the product of charge 
exchange has been sho\m by Hundhausen (196S) (197“^) to be 

unimportant in the problem, the source term involves only losses 

of primary interstellar hydrogen as the result of both photoionization 
and charge exchange. Then the source term is 


^ = - («p+ 


“ “ (tf + or n) f 

' pe 2 o ] 


N 


(2.3-2) 


with rate constants defined in sections 1.4a and 1.4b. 

We note that equation (2) depends on the proton number density n 
of the plasma flow problem. Since the hydrodynamic equations of the 
latter problem depend on f^ through their source terms, the two problems 
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are coupled. But the solution of the plasma flow problem is at the 
forefront of computational capability even to the extent the distri- 
bution of hydrogen is known. Thus. hope of solving the combined problem 
hinges on the separation of the two. 

We have seen in section 1,4c that the source terms represent only 
a first order perturbation on the plasma flow problem. Hence the ap- 
proach \>je adopt (at least for a first iteration) is to assume the unper 

turbed, extended Parker solution for n , That solution we recall as 
2 


n = 


6 

n —S’ in the supersonic region 

G 


4 n — 7 T rn the subsonic region, 
e 2 
r 

s 


(2.3-3) 


Finally because the peculiar velocity is an independent variable 
we have 


div V = 0 


(2.3-4) 


Thus the hydrogen transport is governed by 




(2.3-5) 


with the approximations (2) and (3). - 

As the result of the assumptions, the trajectories of hydrogen 

t 

into the solar wind are straight lines and equation (5) may b^. inter- 


preted 


d f S 
N ^ 

d^ V 


(2.3-6) 



Ill 


where v = -^ and 6l is an infinitesmal element of length along a 
dt 

trajectory. 

In the light of equations (2) and (3), calling 
dropping the subscript N we have 

M = ^ . . (2.3-7) 

f V 

Formally the solution of equation (7) is 


f = f exp “ f ~ 
o J V 

I 


(2.3-8) 


where f is the unattenuated distribution function which in the rest 
o 

frame of the gas is Maxwellian. 

In order to carry out the integrations implied by equation (8) 


w-.*-a-iUC*.w>uOno V 4 aIIU 


seme gecmetrical con 


siderations. As in the case of the extended Parker solution, we take 


the positive x axis of a heliocentric cartesian reference frame 
(X , y , z) to have the same direction and sense as the velocity vec- 
tor V of the sun with respect to the local interstellar gas. With- 
.,.0 

out loss of generality we consider a point in the x , y plane, the 
point being described by the conventional polar coordinates r , 9 , 
Some geometric relations governing linear trajectories of primary inter 


stellar hydrogen atoms passing through such a point follow. 

The peculiar velocities of the hydrogen atoms with respect to their 
mean velocity - V is the vector field c represented in conventional 
spherical polar coordinates as c , I] , | based on the positive x 
direction. The velocity v of a hydrogen atom relative to the sun is 



then given by 


(2,3^9) 



(2.3-9) 


= (c cos T] - V ) + c sin t) cos 5 + c sin T] sin 5 

~ i 


If p is the elevation angle the trajectory of a hydrogen atom makes 
with the X direction in passing through the point r , 0 then (Figure 2.) 


tan 8 = 


c sin T] 

V - c cos T| 
o 


(2.3-10) 


The trajectory vector r of such an atom is described by the coordinates 


X - r cos 6 + R cos 3 
^ = r sin 0 + R sin 3 % 

z = R sin 3 sin 5 


V ^ ✓ 


where R is the distance of a point on the trajectory from the point 
r , 0 and, in terms of the angle coordinate % of the peculiar velo- 
city of an atom describing the trajectory, the azimuthal trajectory angle 
X is given by 

(2.3-12) 

The distance of the sun from a point on the trajectory is 
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In equation (13) the angle y lies at the apex between the incoming 
trajectory to the point r , 0 and the prolongation -of the line from 
the sun to the point r , 0 Finally, in obtaining the magnitude r 
from equation (11), one arrives at the useful identity 

cos y = cos 9 cos p + sin 0 sin 3 cos ^ (2.3-14) 

As a practical matter we do not compute f(r,0) in the helio- 
centric reference frame since we need only its moments for the source 
terms (section 2.2a) of the hydrodynamic equations. Rather we compute 
the contributions to the moments of f(r,0) from atoms in the distri- 
bution expressed in the rest frame of the hydrogen gas. Then to get 
the moments we sum the contributions over the coordinates c , 1] , % 
of the distribution expressed in the rest frame. 

For atoms in the differential velocity space element 6c,.6T| , 

the contribution 6N to the zeroth moment or number density N 

o o 

at large distances from the sun is given according to hypothesis by the 
Maxwellian distribution as 


6N = N c^/c^ exp ( - c^/c^) sin T] 6c 6T] 6^ . (2.3-15) 

o O s s 


In the heliocentric reference frame such atoms have the velocity given 
by equation (9) and the associated speed 


2 2 1/2 
V = (c - 2c V cos 'll + V ) ' 
o o 


(2.3-16) 


along the trajectory whose parameters in the heliocentric frame are 
given by equations (10), (11), (12), (13), and (14). Then in the 
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vicinity of the sun where there has been some perturbing effect due to 
pho toionlzation and (if the trajectory has intersected the solar wind 
cavity, section 1.4) charge exchange, the contribution of the remaining 
atoms 6N to N will be given, from equation (8), by 

J 

6N = 6N exp - P — d^' (2.3-17) 

o V 

I 


Evidently from the discussion above the parameters of the line integral 

over the attenuating path -t s.re all functions of c , T| , ^ as is 

6N and, hence, 6N . Substituting for a ^ (X + Of from equation 
o pc 

(2), we have 

2 - ^ 

(X a » ' Of n 

S ( 2 . 3 - 18 ) 

I I I 

with n given by equation (3) . 

To carry out the integrations explicitely we choose as the dummy 
variable of integration the distance R of the trajectory point from 
the end point r , 0 . Then, since photoionization acts everyt^here, we 
have for the first integral on the right hand side, from equation (13), 


- S 


V‘e 

^ r^ 


2 

Of r 
pe e 


0 

I r 


dR 


2 ^ - 2 
R +2rR cos y +r 


2 

-cr r y 
pe e * 

vr sin y 


( 2 . 3 - 19 ) 


In the integrated expression y in the numerator may be seen to be the 
angle subtended at the sun by the entire trajectory to the end point 
r , 0 ; and in the denominator r sin y is the closest approach of 

the trajectory (or its prolongation beyond r , 0) to the sun. 
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The second integral in equation (13) in general may* contain integra- 
tions over segments of the trajectory intersecting both the subsonic and 
supersonic domains of the solar wind. In the case of the subsonic region 

(II), the solar wind proton number density n is constant, from equation 

* \ 

(3); and the contribution to the integral involves simply the total geo- 

^ i 

metric path length PL II in the region 


O' n 



III 


4a n r 

— (PLII) . 
vr 

s 


(2.3-20) 


In the supersonic region (I) , the solar wind proton density varies as 
-2 

r as has been the case treated for photoionization and the contribution 
to the integral (18) is 


_ 2 

an a n r 

r — di' = — ^ - ^ 

J V vr sxn -y 


(^I) 


(2.3-21) 


where M.I is the angle subtended at the sun by the segment of the 
trajectory intersecting the supersonic region. The angle is given by 


cos ^ + Y 

r * 

s 


H 

2 


MI == 


2 cos 


-1 r sin Y 
r 


r < r 


I - 'J 


(2.3-22) 


To summarize the integrations we have 


5N(c,71,Q = 6N^(c, 71,§) exp - 


2 

a r Y 


vr sin Y 
2 


4a n^r a n r 

+ - (Pbll) + ° ® ® 

V r 
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vr sin Y 


(MI) 


] 




3-23) 
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With equation (23), the moments of f(r,0) required for the mass, 
momentum, and energy soutce terms in the hydrodynamic equations, see 
equations (2,2a-26a, 27b, 28b), are 


N = J f d^v = JJJ" 6N(c,T1,5) 

J vfd^v = J*// v(c,Tl,§) 6N(c,Tl,|) 

/ v^f d^ v^(c,Tl,Q 8N(C,T1,5) 


(2.3-24) 

(2.3-25) 

(2.3-26) 


\ 



118 


Chapter III. Numerical Results - 
Interplanetary Neutral Hydrogen and the Speed 
of the Interstellar Wind 

In this chapter we describe and present some particular results 
from a numerical algorithm developed to compute the moments of the 
distribution of primary neutral hydrogen and the associated source terms 
in the hydrodynamic equations. Specifically, from computed neutral hy- 
drogen number density profiles on the upstream and downstream axes of 
flow we further compute the intensity of backseat tered solar Lyman-of 
radiation expected to be observed at 1 a.u. The computations have been 
carried out for a varied set of conditions with respect to number density, 
temperature, and velocity of the local interstellar gas. Finally the 
theoretical predictions represented by our computations are compared 
with measurements of Lyman-a radiation from earth satellites. Given 
current estimates of the temperature of the interstellar gas, the com- 
parison suggests a meaningful estimate of the speed of the interstellar 
wind . 

1. Algorithm for Computing Moments of the Hydrogen Distribution 
Function. ^ 

The algorithm described in this section implements the approach 
described in section 2.3, We retain here the geometry and notation 
previously defined. 

The algorithm makes use of the fact that hydrogen atoms whose tra- 
jectories pass through the point r , 0 and that have the same peculiar 
velocity coordinate all have trajectories lying in the same plane | = 
f - TT that contains a line through r , 0 parallel to the x axis. 

The given line is viewed as a polar axis of rotation and the azimuthal 
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angles I are measured with the y axis in the y , z plane. For 
evident reason the planes ^ constant for 5 constant are termed 
trajectory planes for the point r , 0 • Within a trajectory plane 
a given trajectory line is described by the angle ^ , Figure 2. 

From equations (2,3-10 and 16) both the angle of the trajectory 
line and the speed of an atom along the line depend jointly on the pecu- 
liar velocity coordinates c and 1j , i,e, a multiplicity of pairs 


c , 1) find trajectories along the same trajectory line. For a given 
trajectory characterized by c and Tj , the attenuation path length 
the argument of the exponential attenuation factor in equation (2,3-23) 
— is seen to be the product of the reciprocal of the speed v and a 
function composed solely of geometric factors depending only on r , 0 , 
P , and ^ , the geometry of the plasma flowfield being given. The 
geometric function which we term the velocity path length contains the 


essence of the prcblem and is costly to compute. Accordingly, since the 


velocity path length is the same for the multiplicity of trajectories 
that can fall along a given trajectory line, for each point r , 0 and 


computational plane § the velocity path length is pretabulated at 
appropriately fine P intervals for later look-up and interpolation as 


a function of P(c,'I|) during integration. 


The computation of the moments of the distribution function is per- 
formed by integrating equations (2.3-24,25, and 26) with equation (2.3-23), 
in order over Tl , c , 5 ( 1 ) . These independent variables are of course 

quantized and the integrations are performed using Simpson* s rule with 

/ 

/ 

the output of the present level of integration being tabulated! 'as argu- 


ments for the next. 
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For each § plane, in the process of integrating the contribu- 
tions to the worttcnts of f the speed v and the velocity v are com- . 
puted for each coordinate pair 1] and c from equations (2,3-9 and 16). 
The angle 3 is computed from equation (2.3-10) and the velocity path 
length obtained from the interpolation table as suggested above. The 
attenuation path length for the trajectory is then simply the velocity 
path length divided by the speed v and the associated e^xponential atten- 
uation factor is computed. 

The initial unattenuated number density element 6N which appears 

o 

as the coefficient of the attenuation factor in equation (2,3-23) is 
given by equation (2.3-15) where it is seen to involve only the product 
of a function of c by a function of Tj , both functions independent of 
^ , For the evident economy of computation, these homogeneous functions 
have also been precomputed and tabulated for the quantized values of 
Lheij. aiguiueiiLs c auu 1] used iu Llie iiiLegrat ioii. 

In computing the interpolation table for velocity path length, the 
figure of the heliopause boundary is approximated by an ellipsoid - 
cylinder, Figure 3. The sun is at the focus of the ellipsoid which fits 
the extended Parker solution exactly at the stagnation point and again 
in a plane that includes the sun and is normal to the axis of flow. 
Downstream of the sun, the cylinder terminates the ellipsoid smoothly at 
the semiminor axis providing another fit to the extended Parker solution 
asymptotically. 

Making available the analytical geometric power of quadratic sur- 
faces, the ellipsoid cylinder approximation permits closed analytic 
expressions for a number of geometric parameters including the length 
of the trajectory segment in the subsonic solar wind region, equation 




Figure 3, Ellipsoid cylinder approximation to the figure of the heliopause boundary. 


121 


122 


(2,3*20). j[n toto the approach results in closed analytic expressions 
for the velocity path length and obviates otherwise prohibitively costly 
numerical integration along trajectories, such integration involving 
searches for intersections with the heliopause boundary surface. 

* 5 * 

i 

2. Axial Profiles of Hydrogen Number Density 

With reference to the algorithm of the previous section, the results 
presented below were computed with the angles T| and § each quantized 
in 3^, increments over the interval 0° to 180^ (converted to radians). 

The speeds c were quantized in 0.1 c^ increments over the interval 
0 to 3.5 c where c is the most probable thermal speed in the Maxwellian 
distribution. The velocity path length table had 129 entries spaced equal- 
ly over the interval 0 to tt . 

Figures 4 and 5 show profiles of Hydrogen number density computed on 

the upstream and downstream axes of flow in the solar wind. In both fig- 

_ 3 A o 

ures the results shown are for N ~ 0.1 cm and T = 10 K, Figure 4 

o o 

is for = 16.6 km/sec and Figure 5, for 8kjn/s. The former speed is 

sonic with respect to the ionized component of the interstellar wind. 

The latter subsonic speed is typical of the estimate we determine in the 

next section from comparing additional computed results with experimental 
\ ■ . ■ 
measurements of Rayleigh backscatter of solar Lyman-a radiation. Also 

plotted for comparison in Figures 4 and 5 are curves derived from an 

analytical formula for the case of zero temperature hydrogen penetrating 

along the upstream axis of flow. 

The flowfield described by the extended Parker solution (section 1,2) 

/ 

and to which Figure 4 refers has the shock at 123 a.,u. and the stagnation 
point at 182 a.u. Similarly for Figure 5 the respective numbers are 



142 a.u, and 282 a.u. 




Figure 4. 


Hydrogen number density profiles 
axes of flow for N «= 0,1 cm“^ 


computed on the upstream. and downstream 

V = 16, 6 km/so 
o 


9 


123 



Figure 5, Hydrogen number density profiles computed on the upstream and 

downstream axes of flow for N « 0,1 cm“3^ V = 8 km/s* 

0 0 
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Looking first at Figure 4 and comparing the number density profiles 
for hot and cold hydrogen on the upstream axis, we find close quantitative 
as well as qualitative agreement for distances beyond 10 a.u. from the 
sun. Within 10 a.u., however, the predicted surviving density of hot 
hydrogen is systematically higher than for cold. 

The trend x^ithin 10 a.u. is easily understood particularly when one 
notes that the bulk speed and most probable thermal speed are comparable 
and where the trend becomes pronounced the surviving number density has 
fallen to roughly half the far field value. The interpretation is based 
on the fact that the peculiar velocity components for the fraction of 
hot hydrogen in the forward T| cone (T| < of the distribution tend 
to reduce the speed v of penetration along the trajectory and thus in- 
crease the attenuation while the converse is true in the rearward 1] 
cone. Thus after much of the slower hot hydrogen in the forward cone 
has been lost to at cenuauion the surviving hot hydrugeu is perictrating 
relatively faster than the surviving cold hydrogen and thence penetrating 
deeper into the solar wind before suffering photoionization or charge 
exchange. The same trend is apparent in Figure 5 but there because the 
peculiar velocities are relatively higher than the bulk speed the effect 
is considerably enhanced, ^ 

In Figure 6 we observe the steep drop in hydrogen number density in 
passing through the relatively proton dense subsonic solar wind region 
on the upstream axis. Between the two curves for the upstream axis, 
the drop in number density is greater through the thicker subsonic region 
associated with the slower bulk velocity = 8 km/s . j 

On the downstream axis the existence of any neutral hydrogen is the 
result of the transverse peculiar velocity components associated with 
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temperature. The higher the thermal speed relative to the bulk speed 
the steeper the angles p for trajectories and the 'closer to the sun 
can hydrogen atoms fill in the cavity on the downstream axis. Farther 
from the sun on the downstream axis the mean penetration distance through 
the subsonic region and the mean penetration speed govern the surviving 
number density jointly. The former decreases while the latter increases 
with leading to increased penetration of the solar wind and in 

Figure 6 the result that the tv;o number density profiles for the rear- 
ward axis eventually cross. The existence of maxima and subsequent de- 
cline in both profiles with distance on the rearward axis is associated 
with the evasion of the denser solar wind regions near the sun on the one 
hand and the progressive lengthening of attenuation paths through the 
subsonic solar wind region on the other. 

Finally in Figure 6 we observe that for hydrogen at the same tem- 
peracure a difference in assumed bulk speed has a considerably greater 
impact on predicted surviving number density on the foreward axis than 
on the rearward. 


3. Predicted Intensity of Solar Lyman-o' Backscatter from Hydrogen 

When the intensity of solar Lyman-cv backscatter at 1 a.u. is mea- 

^ 10 ^ -2 -1 -1 
sured in Rayleighs (1 Rayleigh ^ photons cm sec ster ) the 

derivative of the intensity and the intensity (section 1.4d) are re- 
presented on the axis of flow by . 
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(3.3-1) 
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(3.3-2) 
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where H is the hydrogen number density, the single atom scattering 

-3 -1 13 

rate at 1 a.u. has the value 2.28 X 10 s , is 1.5 X 10 cm , 

and r is measured in a.u. For the number density profiles shown in 
Figure 5 for hot hydrogen, the profiles of the derivative of the inten- 
sity are plotted in Figure 7. We note that the largest contributions 
to the predicted intensity come from distances closer to 1 a.u. on the 
upstream axis than on the downstream axis and that the intensity integrals 
in both directions find effective convergence in the vicinity of 100 or 


200 a.u. . 

Table 7 contains relevant parameters and results from a varied set 

of computer runs made with the model. In the table U and D stand 
respectively for the predicted intensities at 1 a.u. on the upstream and 
downstream axes of flow and U/D stands for the ratio of these intensities 
Inspection of the table reveals the following qualitative trends 
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(3.3-3) 


The trends evident in the last column in relations (3) make the 
ratio U/D an attractive estimator of the interstellar wind velocity 
from satellite measurements of Lyman-a backscattcr. The ratio of 
'^asurements^is independent of error in the absolute calibration of the 
measuring instruments. And the ratio of predicted intensities are 
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Table 7. Synopsis of computer runs with model predicting 
solar Lyman~cy backscattcr on upstream and down- 
stream axes of flow 


T 

N 

V 

U 

D 

U/D 

INT. 

R 

X 

o 

(°K) 

o 

(era 

■ o 

(km/s) 

(Rayleigh) 

LIM. 

(a.u.) 

s 

(a.u.) 

o 

(a.u.) 

10^ 

0.2 

16.6 

1510 

339 

4.45 

100 

98 

129 

ft 

II 

6 

788 

374 

2.11 

II 

124 

256 

11 

II 

0.6 

312 

205 

1.52 

It 

129 

779 

10^ 

It 

6 

1604 

1187 

1.35 

II 

55 

161 

4 

4x10 

II 
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I Figure 8* Plot of velocity of the Interstellar wind vs* computed upstream to downstream 
! Lyman-a backscatter intensity ratio for parameter temperature* 
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Insensitive to the estimate of the density of the interstellar gas. 

Thus the estimator of the interstellar wind velocity from the U/D .• 
ratio is sensitive mostly to the estimate of the temperature of the 
Interstellar gas. 

From the data of Table 7 the predicted ratios U/D are plotted 

versus interstellar wind speed in Fi^re 8. The points plotted with 

an 3^ are for integrations cut off at 100 a.u. and with a + at 

200 a. u. Where two of more points are for the same temperature, curves 

of constant temperature are drawn. Evidently substantial portions of 

4 o 3 o 

the constant temperature curves for T =10 and T = 5 X 10 k 

o o 

approximate straight lines which intercept the point (1,0). Theoreti- 
cally the point (1,0) would be occupied in the instance of the spheri- 
.cally symmetric flowfield that would obtain were the sun at rest in the 
interstellar gas. Straight lines to indicate both the trend of the con- 
stant temperature curves and the spatial organization of the data with 
temperature have been drawn through the points - even single ones - for 
all temperatures represented in the set of computations. 

Because the angles that trajectory lines make for different 

and C are geometric invariants for constant ratio V /c (or equi- 
s o s ' ^ 

valently A/ T~ ), it is of interest to replot the data of Table 7 

and Figure 8 as the ratios V /c vs U/D . Figure 9 shows the data so 

os 

plotted. With the exception of one point the data so organized describe 
a universal curve for which V^/^/ can be regarded as the independent 
variable. 

For the points which fall on the curve, the dashed line indicates 

the curve is highly linear for V /c <1 , The point for V =0.6 km/s 

0 8 . . o 

4 o 

T = 10 K falls off the curve but in the region of the graph 
o 

« 1, U/D 1) where the curve is linear. The point for 




I 


Figure 9* Universal plot of ratio of interstellar wind velocity to most probable thermal 

speed vs, ratio of computed upstream to downstream Lyman-a backscatter intensities. 
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= 16,6 kra/s, = 10 K falls on the curve but out of the region of 
linearity (V /c ^ !)• The evidence is that the model predicts results 

O S • 

that describe a curve that is linear and will occupy the point (1,0) only 

in the limit T , V large and finite but that such a curve is not 

o o 

described for T finite and V small approaching the limit V “* 0 . 
As can be seen in the model for the extended Parker solution, the condi- 
tions securing the domain of linearity above have to do physically with 
restricting the size and relative importance of the subsonic region of 
the solar wind in favor of the supersonic region. For small the 

aspherical subsonic region becomes enormous (Table 7). 

4. Estimate of the Mach Number and Speed of the Interstellar Wind 
Table 8 contains intensity ratios I /l , obtained from 
satellite measurements reported by Bertaux and Blamont (1971, 1972) 
and Thomas and Krassa (1971). (See Thomas, 1972.) The authors cited 
have argued the correspondence of their observations to the physical 
model to which we have given the mathematical representation and from 
which computed the predicted intensity ratios u/D above. In particular 
and in conformity with the predictions .of our model, they have inferred 
that the maximum and minimum intensities should be observed respectively 
when the observational lines of sight are along the upstream and down- 
stream axes of flow. 

Table 8 - Maximum and minimum iotensity^ 

of solar Lyman-a backscatter measured 
by earth satellite at 1 a.u. 


Maximum 

Minimum 

Ratio 

Instrument 

580 

240 

2.42 

LASP 

570 

250 

2.28 

LASP 

460 

180 

2.55 

Paris 

525 

200 

2.62 

Paris 


^data extracted from Figure 3 of Thomas (1972) 
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We note that the measured intensity ratios fall in the linear portion of 
our universal curve which accordingly we may represent by the expression 

V /c = (U/D - 1) . (3.4-1) 

Since the speeds of sound in either the .neutral or ionized com- 
ponents of the interstellar gas and the most probable thermal speed in 
the neutral gas are all proportional to , the left hand side of 

equation (3,4-1) is proportional to the Mach number in the continuum 
flow of the ionized fraction. Thus from equation (3,4-1) with the 

inferred correspondence between I /l . and U/D and from measurements 

max min 

of I /l , we derive the direct estimator of the Mach number of the 
max min 

sun in the' interstellar wind 

(3.4-2) ' 

The mean of the intensity ratios in Table 8 is 2.47 for which 

equation (2) gives the Mach number 0.475. Rounding off meaningless 

digits, we take 0.5 as our estimate for the Mach number of the sun in 

the local interstellar medium. 

4 o 

Further taking 10 K for the estimate of T^ we have 12.8 km/s 

for c . Then from equation (3,4-1) with U/D replaced by the estimate 
_ 8 

2.47 for I /l . we get 7,84 km/s for V , Thus we take 8 Idn/s as 
max min ^ o 

our estimate of the speed of the sun with respect to the local interstel- 
lar gas. 

This number differs substantially but not radically from the earlier 
estimate of 6 km/s found for the same assumed temperature (10^ ^K) by 
Thomas, 1972 with a less detailed model, Thomas * model featured an inter* 
stellar wind having a single constant speed, isotropically distributed 
thermal component traversing a constant speed supersonic solar wind of 



infinite extent 
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Chapter IV - Sunnnary and Discussion 
The termination of the solar wind has been described from the 
theoretical point of view as the problem of finding a model for the 
extended solar corona such that admissible boundary conditions match 

conditions in the perturbed local interstellair medium 

" i 

As a prelude to formulating a mathematical model, certain theoretical 

foundations of our understanding of the problem have been reexamined in 

Chapter I« in the light of most recent findings concerning the states of 

the solar wind and the local interstellar medium. In particular, the 

findings that have had the greatest impact on the present exposition and 

have served in considerable measure to distinguish it from earlier work are 

A o 

that the local interstellar medium has a high temperature ^ 10 K, is par- 
• tially ionized «,10%, and blows past the sun with subsonic velocity ^ Mach 0.5. 

Given these findings, investigation of Parker *s early (1961) analyti- 
cal model for the termination of the solar wind suggested that a simple 
reinterpretation or extent ion of the model provides a useful approximate 
description of the adjacent flows of the solar wind and the ionized inter- 
stellar wind. 

An investigation of solar wind temperature along lines of the one 
dimensional compressible flow equation and heating associated with poly- 
trope solutions brought to light additional evidence that either or both 
the classical coefficient of thermal conductivity is too high — perhaps 
by as much as a factor of five — and the associated classical inhibition 
of heat conduction by the embedded transverse magnetic field is overstated 
for the bulk of low temperature («.10^ °K) electrons comprising the quiet 
solar wind in the vicinity of earth. 

The scale of effects of photoionization and charge exchange are investi- 
gated both for protons and neutral hydrogen atoms in the three regimes of 
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plasma flow — supersonic solar wind, subsonic solar wind, and ionized 
interstellar wind. In the analysis the properties of the plasma flow 
regimes were taken from the extended Parker solution. The results 
suggest that in the ionized interstellar wind region photoionization and 
charge exchange may be ignored in the solution of the hydrogen transport 

. j 

problem but that charge exchange does contribute a significant drag on 
the flow of the ionized component of the Interstellar wind around the 
heli'opause. In the subsonic solar wind region photoionization is not 
important; charge exchange, however, while having only small effect on 
the transport of hydrogen, acts significantly to cool the flow of the solar 
wind. In the supersonic solar wind both photoionization and charge 
exchange have an impact on the transport of hydrogen and on the flow of 
the solar wind. 

In the latter region a linear perturbation solution has been found 
that adequately describes the effects of both processes on the solar wind 
flow. The solution exhibits the point that as the result of photoioniza- 
tlon of hydrogen the addition of mass and associated decrease of velocity 
of the supersonic solar wind has no first order effect on the momentum 
flux and thus the location of the shock surface. The effect of charge 
exchange was found to decrease by 14% the estimate of the shock radius 
as given by the extended Parker solution. 

The preliminary investigation concluded with the finding that within 
experimental uncertainties the opposing forces of solar gravitation and 
radiation pressure from solar L37man-(y scattering are in balance for hydro- 
gen. Then, in view also of the large mean free paths for collisions of 
all kinds, the assumption was made that the trajectories of interstellar 
hydrogen atoms are straight lines into the solar system. 
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In Chapter II. a mathematical model was formulated to describe the 
transport of Interstellar hydrogen into the solar wiftd and the adjacent 
flow of the solar wind with the ionized component of the interstellar 
wind. The global formulation begins with the Boltzmann's equations for 

't; . 

both neutral and plasma species, ; 

yf 

ultimately the flow of the plasma species is described through a 
model on the level of one-fluid magnetohydrodynamics. The model is 
developed through the plasma transport equations in the manner of and 
with the transport terms given by Braginskii (1965). Source terms for 
the mass, momentum, and energy equations in conservation law form derive 
from the interaction of the neutral component of the interstellar wind 
with the solar wind and ionized interstellar wind plasmas. Except in 
the immediate vicinity of the stagnation point singularity, the frozen 
magnetic field approximation has validity. In the solar wind, quasi- 
normal, strong shock relations apply on a closed surface surrounding 
the sun. Tangential discontinuity relations apply at the free surface 
(heliopause) between the solar wind and ionized interstellar wind flows. 

We note that the latter relations demand that the normal component 
of magnetic field vanish at the heliopause. Except in the unlikely case 
of .^purely field aligned flow, the vanishing of the normal field component 
will lead to an imbalance of magnetic field strength around the perimeter 
of the intersection of the heliopause with any plane normal to the axis 

of flow. The associated imbalance in the Lorentz forces of pressure nor- 

% 

mal to the perimeter and tension along the perimeter must then progressive- 
ly push in the perimeter where the field is concentrated and permit the 
perimeter to squish out where the field is weakened, the motion taking place 
in time as the plane and the perimeter are convected downstream with the flow. 



Thus the shape of the heliopause downstream of the sun is not expe.cted to 

be axisymmetric as described by the extended Parker solution but, rather, 

to be “delta winged”, tapering from more or less circular cross sections 

near the sun to thin ellipse-like sections far downstream. Of course, the 

asphericity of the near interplanetary magnetic field embedded in the solar 

.i 

wind must also contribute to the departure of the heliopause from axis 5 nn- 
metry. To describe in some quantitative fashion the scale of the depar- 
ture of the heliopause from axisymraetry must be regarded as one of the 
outstanding theoretical problems associated with the termination of the 
solar wind. 

In the field free, single collision approximation adopted for the 
transport of neutral hydrogen, the Boltzmann equation reduces to a con- 
tinuity equation for the distribution function. Since the fast secondary 
hydrogen products of charge exchange have been shown by Holzer (1972) to 
have negligible importance, the continuity equation has only source terms 
for removal of hydrogen through photoionization and charge exchange. But ► 
the source term for the neutrals to suffer charge exchange depends on the 
density of ions and the source terms of the hydrodynamic equations for the 
plasma flow depend on the density of neutrals. Thus the two formulations 

are coupled. To resolve the impasse at least for the first iteration of 

\ . 

solution, the extended Parker solution was assumed for the plasma and the 
hydrogen transport problem solved in the approximation afforded. 

To the latter end, an efficient semianalytical numerical algorithm 
(Chapter III.) has been developed to carry out the solution on the needed 
level of moments of the hydrogen distribution function. The algorithm 
depends on approximating the figure of the heliopause by a prolate ellip- 
soid — cylinder which is fit to the heliopause, first, at the stagnation 
point; secondly, in a plane that includes the sun and is normal to the axis 
of flow; and, .finally, assymptotically far downstream. 
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The algorithm has been used to compute profiles of hydrogen number 
density on the upstream and downstream axes of flow.' Using the single 
scattering model developed at the end of Chapter I., we have computed 
intensities of backscattered solar Lyman-a radiation predicted to be 
observed at 1 a.u. as the result of scatter from hydrogen distributed 
according to the previously computed axial profiles. The results of 
computations for a varied set of parameters of the interstellar wind has 
been found over a substantial range to exhibit a strong linear dependence 
between predicted ratios of upstream to downstream backscatter intensity 
and the ratios of the independent variables bulk velocity and most prob- 
able thermal speed of the interstellar hydrogen. 

Bertaux and Blamont (1972) and Thomas (1972) have provided the inter- 
pretation that the maximum and minimum backscatter intensities found in 
satellite observations occur repectively when the lines of sight are along 
the upstream and downstream axes of flow. Since the ratios of maximum to 
minimum intensity observed by satellite fall in the linear .range of our 
predictive model for upstream to downstream intensity ratios and given 
the inferred correspondence between these ratios, we inferred that obser- 
vation of maximum to minimum backscatter intensity provides a sensitive 
measure of the ratio of the speed of the sun through the local interstellar 
gas to the most probable thermal speed of the gas. Further, since the 
most probable thermal speed and the speed of sound are in constant ratio, 
the observation of Lyman-o backscatter provides a sensitive measure of the 
Mach number of the sun through the interstellar gas. Based on the ratios 
of observed maximum to minimum backscatter intensity reported by Bertaux 
and Blamont (1971) and Thomas and Krassa (1971) the Mach number that we 
found is 0.5 with respect to the ionized component of the interstellar wind. 
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Given the estimate 10 K for the temperature of the interstellar medium, 
the corresponding estimate for the speed of the interstellar wind is 8 km/s. 
We note that with respect to Mach number, consistency has been found 
between the predictions of the model and the assumption that the flow is 
subsonic as required by the extended Parker solution. While being fully 
cognizant it is only an approximation, the latter solution we have repre- 
sented as being a useful approximation for the analytical purposes we have 
made of it. In particular since the backscatter intensity integral finds 
effective convergence well within the distance to the stagnation point 
tod most of the scattering comes from within the supersonic region, we do 
not expect anticipated departures of the solution from that used here to 
greatly affect the estimate of the Mach number, 

' However, as Holzer (1972) and Axford (1972) have discussed, the 
strengthening of the magnetic field lines by stretching and curvature in 
both the subsonic region of the solar wind and in the interstellar wind 
and by cooling of the subsonic solar wind through charge exchange introduces 
additional compression in the region beyond that predicted by the Bernoulli 
eqnatloh; To assess the scale of these additional, highly nonlinear, 
three dimensional effects on the dimensions of the flowfleld really 
requires a full numerical solution of a model on the order of that 
developed here. Given present uncertainty in the measurement of astro- 
physical quantities, such a solution may not provide measurably greater 
Information than we currently possess about the absolute state of the 
system. And obtaining such a solution will be a formidable and costly 
task, but having a solution for even one set of nominal conditions — 
hopefully with a transverse interstellar magnetic, field — would consti- 
tute a most fascinating contribution in structural detail. 
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